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ABSTRACT 


A two-level, quasi-geostrophic atmospheric model on a 
mid-latitude ß-plane with a single wave in the x-direction 
was coupled with a three level primitive equation ocean 
model through the exchange of sensible heat and momentum 
to investigate the effect of the ocean circulation on tran- 
sients in the atmospheric circulation and vice versa. A 
long-term (79 years) simulation of the atmosphere-ocean 
system was carried out. The experiment was conducted in 
three phases that consisted of integrating the combined, 
the ocean-only and the atmosphere-only model equations. 

The combined model produced realistic patterns of both 
atmospheric and oceanic circulations. Energy studies of 
the atmosphere were completed with the results comparing 
favorably with the existing knowledge of the energy balance 
in the real atmosphere. 

Data were retained at selected points on a daily basis 
and spectrally analyzed for a simulated six year period. 
Transients with periods on the order of one and two months 
appeared in the spectra of both atmospheric and oceanic 
variables from the combined model. These transients were 
not observed in the ocean-only or atmosphere-only phases of 
the experiment, and are attributed to the interaction occur- 


ring in the combined model. 
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127 INTRODUCTION 


The use of numerical models to simulate physical con- 
cepts that would otherwise be restricted to limited observa- 
tional studies has advanced the knowledge of both the 
atmosphere and ocean and their large scale interaction a 
great deal over the last several years. But that knowledge 
is far from complete and it is the intent of this work to 
reveal more information on the long term, large scale inter- 
actions and the frequencies of the response between the 
atmosphere and ocean. The primary goal of this study was 
to identify the effects of the ocean circulation on the 
transient (time-dependent) atmospheric circulation, and 
similarly the effect of the atmospheric circulation on the 
ocean. 

The development of a large scale air-sea interaction 
model that would remain computationally stable throughout 
a long-term numerical integration of the dynamical equa- 
tions of motion was required to achieve this goal. It was 
desired that the method by which the integrations proceeded 
be one which allowed the synchronization of the oceanic and 
atmospheric time steps so that a "one-to-one" relationship 
existed between the two media. Many complex numerical 
models of both the ocean and the atmosphere have been 
developed in the past. However, those models all require 


large amounts of computer time because of the very short 





time steps required to satisfy the linear computational sta- 
bility requirement. The model developed in this work had 
fewer degrees of freedom in the atmosphere and this resulted 
in faster integration of the atmospheric equations. It also 
utilized certain approximations and assumptions that permit- 
ted the use of longer time increments in the integrations. 
The necessity of employing a technique to simulate time 
synchronization was thereby avoided. 

Many numerical general circulation models have been 
developed in the last two decades. Phillips (1956) made a 
long period forecast with a two-level, quasi-geostrophic 
model, starting with an atmosphere at rest. Diabatic ef- 
fects were included in the equations as a linear function 
of latitude, and such quantities as the intensity of the 
heating, the value of the vertical stability and the type of 
frictional dissipation, were empirically specified. The 
flow patterns which developed were quite realistic. With 
the A of this model, Phillips demonstrated the 
feasibility of numerical simulation of the atmospheric 
general circulation with a simple model. 

Mintz (1958) proposed numerical general circulation 
experiments that specified a heating function in which the 
diabatic heating of the air depended on the air-sea tempera- 
ture difference. He showed that this heating function could 
be used to study the circulation of the atmosphere on a 
planet covered by either land or ocean or a combination 


thereof. 





In both of the above works, the partial differential 
equations were formulated in finite differences. Lorenz 
(1960) discussed a method of expressing a dependent variable 
as a series of orthogonal functions. An example of such an 
analysis, in one dimension, is the Fourier series, in which 
the orthogonal functions are sines and cosines. The coeffi- 
cients of these orthogonal functions become the new dependent 
variables and a finite number of variables are retained for 
prediction. Bryan (1959) reported on his work with Lorenz 
in which this technique was applied to the variables in a 
two-layer, quasi-geostrophic model. The variables were ex- 
panded into a few leading terms of a spherical harmonic 
series to conduct a simple general circulation experiment. 

Smagorinsky (1963) reported on an extended period numeri- 
cal integration of a two-level, baroclinic primitive equa- 
tion model which had "land" as its underlying surface. This 
experiment was made for the simulation and study of the 
dynamics of the atmospheric general circulation. The solu- 
tions corresponding to external gravitational oscillations 
were filtered by requiring the vertically integrated diver- 
gence to vanish identically. This was a major departure 
from earlier work in that many hydrodynamic and kinematic 
constraints were removed. 

Mintz (1964) conducted an experiment in climate simula- 
tion by numerically integrating the primitive equations of 
atmospheric motion, with friction and heating terms included, 


over a long period of time. His two-level model domain 





encompassed the entire earth, with the underlying oceans 
and continents having different thermal properties. Leith 
(1965) described a six-level numerical model with which an 
attempt was made to compute the hydrodynamic and thermo- 
dynamic evolution of a moist atmosphere on the whole globe. 
The model took into account such external influences as 
solar heating, evaporation and surface friction. 

Smagorinsky, Manabe and Holloway (1965) reported the 
results of a nine-level general circulation model of the 
atmosphere. The primitive equations were used in that 
Study in SL the atmosphere attained a state of quasi- 
equilibrium. The nine levels of the model were distributed 
to resolve surface boundary layer fluxes as well as radia- 
tive transfer by the prescribed distribution of ozone, 
carbon dioxide and water vapor. This model was an outgrowth 
of Phillips' two-level quasi-geostrophic model and Smagor- 
insky's two-level primitive equation model. 

Arakawa (1966) derived a finite difference Jacobian 
expression for the advection term in the difference vorti- 
city equation for two-dimensional incompressible flow that 
properly represented the interaction between grid points. 
He showed that the use of this differencing technique main- 
tained the integral constraints on quadratic quantities of 
physical importance and that this was sufficient to prevent 
nonlinear computational instability. Before Arakawa's work, 


long-term integrations had been possible only with the use 


10 





of large friction and diffusion terms. Arakawa's scheme 
permitted long-term integrations without the requirement for 
such artificialities. 

Kasahara and Washington (1967) described a model of the 
general circulation of the earth's atmosphere which was 
developed at the National Center for Atmospheric Research 
(NCAR). A distinguishing feature of this model was that 
the vertical coordinate was height rather than pressure, 
though hydrostatic equilibrium was maintained in the system. 
Experiments with a two-level version of this model were con- 
ducted. Long-term (over 100 days) integrations were conduc- 
ted testing various physical parameters and examining their 
affect on the results. 

Numerical models have also been developed for predicting 
the atmospheric conditions in enough detail to provide 
daily forecasts for given areas. Kesel and Winninghoff 
(1972) reported on the U. S. Navy operational five-level 
primitive equation model that has been used successfully 
for operational forecasts since 1970. This hemispheric 
model perhaps typifies the state of sophistication to which 
numerical weather prediction has progressed. 

Early mathematical models of the ocean circulation were 
very simple as were the early atmospheric models. Munk 
(1950) utilized the oceanic mass transport to examine a 
wind-driven ocean which did not include specification of 
the vertical distribution of density or currents. He ad- 


hered to the linear theory in obtaining analytic solutions 


Tel 





to She equations for steady horizontal motions. Fofonoff 
loca )eistudted) frictionless, steady horizontal flow which 
could occur in a homogeneous ocean of constant depth. The 
effects of the thermohaline structure and compressibility 
were neglected everywhere. By considering this idealized 
case he avoided the analytical difficulties presented by 
the nonlinear terms in the hydrodynamical equations. 

Veronis (1963a, 1963b, 1965, 1966) conducted numerous 
numerical experiments using both a grid point method of 
examination and a Fourier method of analysis and made con- 
parisons of the two methods. He generated flow patterns in 
a rectangular, B-plane, barotropic ocean by assuming the 
basin could be divided into an interior region in which the 
dynamical balance was quasi~geostrophic and a boundary region 
where the dynamical balance was governed by inertial forces 
much the same as Munk. He developed both linear and non- 
linear models and experimentally determined optimum values 
of parameters, such as the Rossby number and the bottom 
PeletlOmmparameter needed for the numerical work. 

Bryan (1963) obtained solutions for a time dependent, 
nonlinear, barotropic model of a wind driven ocean by 
numerical integration. He varied parameters such as the 
Wind stress, the Rossby number for the interior flow and 
an effective Reynolds number for the boundary current and 
examined the results. He examined boundary currents and 
the separation of the western boundary current from the 


coast for both linear and nonlinear cases. Bryan and Cox 


172 





(967) investigated the ocean circulation using a baro- 
clinic model to numerically integrate the primitive equa- 
tions in almost complete form for different values of 
Reynolds numbers. The near linear solutions for low Rey- 
nolds number were most accessible to a numerical net with 
coarse resolution and gave valuable information for under- 
standing the more complicated regimes at higher Reynolds 
number. 

Gates (1968) integrated the primitive equations for a 
rectangular homogeneous ocean with a free surface on a $f- 
plane starting from an initial state of rest. He investiga- 
ted the transient Rossby waves and found a well-marked life 
cycle of amplification and decay as they propagated west- 
ward at about l m Sean. 

Bryan (1969) continued his earlier work with Cox by 
adding new features to their five-level model. He treated 
salinity separately from temperature and allowed sea ice 
formation in the new model. This work is discussed in more 
detail later in the section on the combined model. Haney 
(1971) discussed the formulation of the net downward heat 
flux at the ocean surface, and used this formulation in 
experiments with a six-level baroclinic ocean model similar 
to Bryan's with emphasis on the equatorial regions. This 
model is representative of the present state of knowledge 
in ocean modeling. 

In recent years research on the problem of atmosphere- 


ocean interaction has progressed to the point where very 


13 





complex atmospheric numerical models have been combined with 
similarly complex ocean models and integrated for extended 
periods of time. The major work in this area combined a 
multilevel, primitive equation atmospheric model including 
the earth's hydrology reported on by Manabe (1969) with a 
multilevel, primitive equation ocean model (Bryan 1969). 
The goal of their study, quite different from the present 
one, was to determine the effect of the ocean circulation 
on the time average climate of the atmosphere. Their ex- 
periments produced some notable results such as realistic 
heat transports and reasonable distribution of the pole- 
ward transport of energy. In one sense, however, this ex- 
periment was considered incomplete in that the goal of 
calculating a climate equilibrium with the combined system 
was not completely attained. 

That research was carried out in three steps: 

a. A study of the atmospheric model without the effect 
of ocean circulation. The effect of heat transfer by the 
ocean current was neglected. The ocean had zero heat 
capacity and simply acted as an infinite reservoir of 
moisture for the atmosphere. 

b. A study of the oceanic model without any feedback 
from the atmosphere. The distribution of surface tempera- 
ture, wind stress and precipitation over the ocean were 


given and held constant with respect to time. 
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a A study of the joint ocean-atmosphere model in which 
the two systems were allowed to interact fully with each 
other. 

A comparison of the results of the three experiments 
yielded some insight into the role of the ocean in maintain- 
ing the climate and of the atmosphere in maintaining the 
thermal and dynamical structure of the ocean, 

In the Manabe-Bryan model more computing time was needed 
to obtain a one day integration with the atmospheric model 
than with the ocean model. Therefore, attempts were made 
to economize computer time by integrating the atmospheric 
equations for only one year while integrating the oceanic 
equations for 100 years. This resulted in a distorted 
response of the ocean circulation to fluctuations of atmos- 
pheric circulation, since an atmospheric disturbance that 
should have an effect for possibly one week influenced ocean 
conditions for up to two years. For example, it took one 
week for a cyclone wave to traverse the oceanic region. 
Because of the synchronization, the oceanic part of the 
model felt that cyclone for as long as two oceanic model 
years. 

Wetherald and Manabe (1972) reported on the effect of 
seasonally varying the solar radiation in the combined 
model just discussed., Many of the results obtained in the 
new version of the model were different, the most signifi- 
cant of which was a warming trend which occurred when 


seasonal variations of the solar radiation were introduced. 
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It should be pointed out that an exact statistical equili- 
brium was not obtained in any part of the region of this 
integration. 

After an extensive review of the literature which in- 
cluded all of the above works, it was decided that to reach 
a State of quasi-equilibrium with a coupled and time- 
synchronous atmosphere-ocean model the simplest atmosphere- 
ocean system which still contained the essential features 
of the dynamics of each media was required. The atmospheric 
model selected was a version of Phillips’ (1956) two level, 
quasi-geostrophic model modified by implementing Lorenz's 
(1960) idea of replacing the dependent variable with Fourier 
coefficients. This decision was based primarily on the 
speed with which the equations of this model could be 
integrated. The equations governing the atmosphere could 
be integrated for one model year in about 30 minutes on the 
IBM 360/67. A model similar to the one desired had been 
used in an earlier atmosphere-only study at the Naval Post- 
graduate School by Taylor (1970) and was readily available 
for further modification. The details of the atmospheric 
model are described in Section II. 

The ocean model selected was a modified form of Haney's 
(197la) six level baroclinic model. The number of levels 
was reduced to three and the horizontal domain reduced to 
an area in the Northern Hemisphere. The depth was reduced 
and it was further simplified by using the ß-plane approxi- 


mation. While this primitive equation model was fairly 
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complex, it was readily available for modification. De- 
tails of the ocean model are described in Section III. 

The ocean and atmosphere models were joined through 
the exchange of momentum and sensible heat. Unlike the 
model of Manabe (1969a), the atmosphere contained no odie 
ture so there was no explicit exchange of moisture. In- 
stead, the effects of evaporation and precipitation 
processes were simulated by other means. These simulations 
as well as the details of the coupling of the two models 
are described in Section IV. 

In a manner somewhat similar to the experiment of 
Manabe and Bryan, the research was conducted in three phases. 
In phase I, equations of the combined model were integrated 
until a state of quasi-equilibrium was achieved. This 
required 79 years of simulation using more than 100 hours on 
the IBM 360/67 computer. The solution of a complex non- 
linear system of equations by numerical integration over a 
long time period as was conducted in this simulation caused 
some concern over the accuracy of the numerical results. 
Error analysis studies by Fisher (1965) and others have 
Rican that, even for short time integrations, errors in 
amplitude and phase arise in the numerical solution. This 
is not to say that such errors make the results meaningless 
because even if fairly large errors were present at the 
end of a long-term simulation of the atmospheric circulation, 
the gross features of the solution would not be changed 


Significantly. However, the results should be examined more 
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from a qualitative than a quantitative point of view. The 
results in this long-term experiment surely contain numeri- 
cal errors but these results are considered to be satisfac- 
tory for this study of the coupled atmosphere-ocean systen. 

One hundred day ss of the atmospheric surface 
stress and surface temperature were calculated and used as 
driving mechanisms for the ocean in phase II, an ocean-only 
experiment. The time of this experiment corresponded to the 
last six years of phase I. Similarly averaged sea surface 
temperatures were used as a forcing function during phase 
III, an atmosphere-only experiment. This also was initial- 
ized so that the time corresponded to the last six years of 
phase I. Details of the procedures used in all three 
phases are described in Section V. 

The resulting data were examined in two ways. The first 
was a synoptic-type study in which maps of both ocean and 
atmospheric parameters were studied and their interactions 
noted. The combined model produced both atmospheric and 
oceanic features that compared favorably with observed fea- 
tures. Energy balance diagrams were constructed for the 
atmosphere and these compared well with existing knowledge 
of the real atmosphere, These synoptic results are shown 
oe Crom) Vil A: 

The second method of examining results was a spectral 
method. Time series of data recorded over the last six 
years for selected locations in the ocean and atmosphere 


were spectrally analyzed to determine if any short or 
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intermediate term fluctuations such as those studied by 
Bjerknes (1962) and Namias (1969) appear due to the inter- 
actions. The results of this analysis are discussed in 
section VI B of this work. By comparing the spectral 
characteristics of the coupled atmosphere with those of 
the atmosphere with fixed ocean surface temperatures, the 
role of the ocean circulation in producing certain tran- 


sient atmospheric responses may be studied. 
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Pi. LHe ATMOSPHERIC MODEL 


A. DESCRIPTION 

The atmospheric model was based on the selective appli- 
cation of the quasi-geostrophic theory to the equations of 
motion. Since the large-scale motions in extra-tropical 
latitudes are quasi-geostrophic in nature, the domain of 
this model was limited to the middle latitudes. The model 
was constructed on a B-plane with a central latitude of 45 
degrees. Phillips (1963) discussed geostrophic motions 
that were characterized by approximate uniformity of poten- 
tial vorticity. These motions have characteristic wave- 
lengths of about 4000 km. He used an expansion in the 
Rossby number about latitude 45 degrees to show that the 
B-plane approximation is justified for ee eee 
of this type. 

Phillips (1956) developed a two-level quasi-geostrophic 
model that he considered was the simplest model in which 
heating and friction could be incorporated. The basic model 
described in this thesis followed Phillips' work in the 
manner in which the quasi-geostrophic equations were derived. 


The quasi-geostrophic vorticity equation on a f-plane is 


e , 2 av dw 
SE (UY) + IVY + (VÍ) + BSE - ES OS 
a OT ÒT 
"talar Tay? =. 
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where TY othe component Of the stress in the Gd direction, 
[> is the component of the stress in the jj direction 
and other symbols have their usual meteorological meaning. 
The model is restricted to two levels with a vertical 
structure as shown in Figure 1. The stream function Ys and 
Y, are representative of the flow patterns in the upper and 


lower layers of the atmosphere respectively. 





O O 

pe ee —250 
2 500 
ee a — — 750 
4 1000 


Figure 1. Vertical structure of atmospheric model. 


Boundary conditions at the top and bottom are: 


w= 0 at p= 0 


and w= 0 at p= Py (2.2) 
The stress Ty vanishes at p = 0. The stress at anemometer 
level was computed from 

eg = p(C,), wi w (2.3) 


2 





where (Ci), 1S a constant drag coefficient and \V is the wind 
at anemometer level. The stress at level four was obtained 
by assuming that \V was a fraction) S, of the vectoriaily exs 
trapolated wind, W,. Phillips (1956) suggested that for 
Simplicity in stress calculations, the anemometer wind speed 
can be assumed to be 70 per cent of the geostrophic wind 


speed W,» The stress was therefore taken to be 
2 
Be De w, I Wy (2.4) 


where Y = 0.7 


A method for estimating E used by Manabe (1969a) is 


h 2 | 
cn Ik /in > (2.5) 


where k Von Karman's constant, 


height at which the stress is to be evaluated 


and 7 


5 roughness parameter. 


The quantity Cy was assumed to be a constant equal to 4.5 

x D. This is larger than values normally used for this 
coefficient, but since Cr depends on where in the surface 
layer the stress is evaluated, the appropriateness of the 
value can best be examined from the resulting momentum ex- 
change between the atmosphere and surface. Values of para- 


meters at level four were obtained by linear extrapolation 


Peep ee the quantity W; was thereby obtained from 


avy. (2.6) 


=- 2 i 
Vie ao We 7 WV, 


4 
Smagorinsky (1963) derived a formula for calculating 


the internal stress at level two, by assuming that the 
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stress was proportional to the vertical wind shear. A 
slightly modified form of this formula was used in this 


work and can be written 
k gp = 
eo ee N, AA) 


where (pk), is an exchange coefficient probably dependent 
on the local Richardson number. An estimate of this from 


Rossby and Montgomery (1935) and used by Smagorinsky was 


50 gm Se ae Other values of this coefficient are 


=l (Palmen, 1955) and 500 gm One ce Reon 


225 gm ee 
1951). the choice of 50 gm oc for this model was 
based on many preliminary numerical experiments with the 


model. The quantity > is the inverse thickness of the 


2 
750-250 mb layer and it is equal to 0.1266 A If typi- 
cal atmospheric wind values are used, the ratio of T, tod, 


is on the order of ten to one. 

The horizontal domain of the model is shown in Figure 2. 
The east-west extent is 4800 km and the north-south distance 
is 6000 km. The width was chosen as 4800 km to coincide 
with the wavelength of the x-wave that will be discussed 
Shortly. The southern boundary is at 18N and the northern 
boundary is at 72N. The use of the ß-plane over such a 
latitudinal extent may affect the results, but only in a 
quantitative manner. 

Mii delos tatie approximation, the equation of state, 


and the geostrophic condition (V > a) were applied to 
O 


the first law of thermodynamics, a quasi~geostrophic 


2 





72% 


6000 


18%N 
4800 


Figure 2. Horizontal domain of the 
model. Distance in kilometers. 


version of that equation would be written 








a (2Y . y(% O E g 
T Cp? + fk x VW N + p Y aan Q , 022,8) 
O O 
where 2 
on DE ea? 
o = — [77 + C Eom sec eb =. 
P 8 P 
= R/C 
s / p , (2.9) 
Q = heating rate per unit mass 
and I: temperature from the standard atmosphere. 
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Dıezabovesconstant value of o, the static stability, was 
used throughout the experiment. 

Equation (2.1) was applied at levels one and three and 
Equation (2.8) was applied at level two. The following 


equations were obtained: 





ow E 
o 2 we: o 2 
Be (VW) tk x O e a r 
a OT IT 
ov, f Ww 
0 2 o 2 
St (V v3) +ik x VY. vv? v5) + BS — + Kp 
oT OT, 
=e T (oe = Su 3 (2,11) 
and 
pe Y,¿+ Y Y OW : 
2 1 3 1 3 1 8 2 K 
E ao o k ae ) as 27 E FEE LP, Q, (22212) 


where Ap = 500 mb. 


Following Lorenz (1960), Var and vn were defined as follows: 


MS VD 
and 2.139 
Y7 Y 
er 


Then Equations (2.10) and (2.11) were added and divided by 


two to yield: 


n 
EV py) tik x Wy, VO P) tik x Ypa e VOY) 


Vu OT Ot 
os Ap E a an 
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Hama ion m2) was solved for Wo» Equation (2.11) was 
subtracted from Equation (2.10) and the result divided by 
two. The quantity w, was replaced and the following was 


obtained: 








9 2 2 ; ae . 2 
TAA - ų dr + fk x vv V(V u vy + lk x vv V(V Vw? 
+ 8 y E wk Q AE en = 
o ox 2f 2 A p ox dy 72 
OT OT, 
se - wf (2.15) 
2 
2f 

-1 -2 

bere EB = D = 3.69 x 10 6 en: 


To reduce computing requirements, the model was further 
simplified by expressing the dependent variables, Yu and 
Ver as functions limited to one wave in the zonal direction. 
The wave chosen was of the predominate baroclinic scale. 
Its wavelength was 4800 km, the zonal width of the model 
atmosphere. Any atmospheric model that is developed with 
the intention of reasonably portraying air~sea interaction 
must include baroclinic disturbances containing wave nun- 
bers of about five to eight. The baroclinic processes that 
are responsible for the resupply of kinetic energy to the 
atmosphere are reasonably represented by the wave number 
chosen. Phillips (1956) and Smagorinsky (1963) have demon- 
strated in experiments with no forcing functions in the 
east-west direction that much of the disturbance energy is 


isolated in a single predominate wave number. 
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By transforming the system of equations to wave-number 
space in the zonal direction, grid points were eliminated 
in the x direction and the requirement for relaxation 
in order to solve for was removed. This type of approach 
had further computational advantages in that no linear 
truncation errors arose in x-space. Also, the fact that 
there were no nonlinear truncation errors in x-space re- 
duced the possibility of nonlinear instability. Therefore, 
errors caused by aliasing due to x-differencing were elim- 


inated. The fields were defined as follows: 


Yy 


ticos kx + B(y,t)sin kx (2,19) 


and 


Ver ISA CNS cos Kx E D(y‚t)sin kx (2.17) 


where A, B, C, and D are Fourier amplitudes of the disturb- 
ance, and E and F are the zonal means. The quantity k is 
the zonal wave number and a value of CO was 
chosen. 

In the initial development of the atmospheric model the 
heating was simply proportional to the air-sea temperature 
difference. A sea surface temperature, I. was specified 
from climatological data and held constant in time. Later 
in the work a more complete form of the heating was used 
and this will be discussed in a later section of this thesis. 

Ihe mean temperature of the layer one-three is propor- 


tional to Y > By using the geostrophic approximation and 


the equation of state, the following equation was obtained; 
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wear nd 
T, R Ver (2.138) 
where T,' is a departure from the mean temperature at level 
two. To obtain physically reasonable values for the atmos- 


pheric temperature at level two a constant was added to T,'. 





Then 
2f 
T, en Yap 525023 a (2.19) 
By using the standard lapse rate, T, was extrapolated down- 
ward to obtain the temperature at level four, T,» and 
oT, 
T,*T, = 577 (42) (2.20) 


where Az is the standard thickness of the 1000-500 mb 
layer and ƏT /əz is the standard lapse rate. Therefore, 


T, only differs from T, by a constant. 


2 
Since the static stability was held constant, the heat- 


ing at level two was proportional to the air-sea temperature 


difference and was expressed as 


Q = —u 1 a 


A p 1 sea 4) ot 


where Hy is a factor derived by Haney (1971) and is a slow- 


ly varying function of latitude. In this work a value of 
55 cal nase was used and held constant in both latitude 
and time. 

When Equations (2.16) and (2.17) were substituted into 
Equations (2,14) and (2.15) and terms of wave number 2k and 


higher were neglected, equations in the Fourier coefficients 


were obtained. The various terms were separated and 
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coefficients of the cosine terms, sine terms and terms in- 
dependent of x were equated. This yielded six predictive 
equations for A, B, C, D, E, and F. The derivation of these 
equations is contained in Appendix A. 


The equations are: 








mo o JE 9%B , 9 39D 9% 93 
Jen BE Sr a 3 ee > 
0 y y oy d y dy 
2,9E oF 
E. 7 2 ea 2 oC _ JA 
+ > scp d w, | k~ (A-2C) re [|w, | (2 A OS 
2 2 2 3 3 
E a, Ey BE 
y oy y oy dy oy 
2,9E oF 
+ k Cy A + Sy C)] + BAk + 
TE arg oD _ 3B 
Ar Ki 20) + ¿LW lO 3, - 3P1) > (2.23) 
2 2 2 
ALL do to oa y) 
dy Y dy 7 dy 2 
3 3 
_ OF oe 1) To E -9 F A 
y en ee O 
2 
re Bao 8 BP, ¿9 € _ ,2 
A A IN A 
O dy 
E F lg? _ ee 2D oe 
aes OSD AE (2D-B) - ¿yLIW,¿ 1 OQ y 7 Dy) 162 +24) 
9 9D 2 2 JE 2 2 IF 2 2 
ee ey teren 0) 
2 2 3 3 
y y + Gar Bel + 8, cx 
dy dy dy dy 
ukg > g 32D 2 go 
k 2f ÊP Pa R 2) Ap a. k US ) 
= 7 Ik? a cr oD _ 9B 
+ E pc, Iw, |x (2D-B) xy LW, | (255 N IES 
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een N 9D 90 
o = > a - BZ + 07 _ Dara) 
dy dy dy dy dy 
Be 9_ oF _ JE 
+ i pc, ay LI, | C2 > a ; (2026) 
2 2 2 2 2 
pela E A rie zA + 0 = 05] 
dy oy dy dy dy 
2 9 UK an 
+ y 3720-aD1| A 
eo pg 3 F B. 
"A A 
9 — a8 IF 
tw, GS = 257)! (2 
where 
OT. a 
GE eu a Wr 


Devcon OF SOLUTION 
l. Finite Difference Scheme 

Centered time differences were used for all terms 
except those involving friction and heating. Those terms 
were computed at time (t-At). A forward time step was 
utilized to start the marching process. A finite differ- 
ence scheme developed by Matsuno (1966) was used every 24 
time steps. This was a two-step iteration to simulate the 
backward difference method. This scheme was applied in 
this situation to eliminate the computational mode that 
would arise from using the centered time differencing ex- 
elusively and also helps eliminate nonlinear instability. 


A time step of 30 minutes was used. 
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Since a wave number representation was used in the 
ereer on, nna te differencing was not required in that 
direction.: However, centered space differences were taken 
Mehendi rection. A uniform grid distance of 200 km was 
used for Ay. 

2. Boundary Conditions 
Free-slip walls existed at the north and south 


boundaries. The boundary conditions at these walls were 


A, B, C, D=0 at y=0, y W (22 8A 


and 
OE 
rn SI Gane oo at y = 0, y = w, (2.28b) 


Condition (2.28a) was imposed to force the meridional flow 


to equal zero and correspond to the statement that the 


normal geostrophic velocity vanished at the wall. From 
Ju 
the zonally averaged equation of motion, — O at the 


north and south walls. The boundary condition (2.28b) en- 
sures that this was true and was also a requirement for 
proper energy conservation. 

tal Conditions 

Initial conditions of a mean zonal current with a 

superimposed disturbance were specified. The vertical 
shear and hence the thermal wind were near the critical 
value for baroclinic instability. The initial disturbance 
had a wavelength, previously discussed, and a structure 
such that when superimposed on a baroclinic current with no 


horizontal shear it had the maximum growth rate. 


Sali 





Hct tal conditions were: 


on sin -3 a Pig 


E = 1200 m cos E ae ee 


and 


2 =a 


F = 600 cos E) cm sec : (2.29a) 


Tt was further desired that the mean surface wind vanish 


Potty oven can be achieved by requiring that 


ay 
a he, eee 


a Teen... = 
4 ay 77 ayitm 7 “pd => pylE-2F1 = 0 EN 


or 


E 2F, 


4. Time Integration Process 

Different methods for solving the time tendencies 
of the Fourier coefficients were utilized. A Gauss-Jordan 
elimination method was used to solve equations (2,22), 
A ZA ZA EZ 25), and (2.27). Richtmyer (1957) dis- 
cussed the use of such a procedure for solving a set of 
equations that had all the elements of its corresponding 
Matrix vanish except those on the three main diagonals. 
These equations had that property. The method of solution 
for Equation (2.27) was discussed by Phillips (1956) and 
modified to ensure that the equation satisfied the boundary 
conditions. Equation (2.26) was solved by a direct march- 


ing process utilizing the proper boundary conditions. 
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Gr 
ive Hadley Regime 

To examine the 
model atmosphere, 
have produced a single 
circulation. A linear 
with warm temperatures 
The amplitudes 

Zero: (A, 


reduced to: 


An assumption t 


PRELIMINARY EXPERIMENTS 


effect of the disturbances on the 


an experiment was conducted which should 


Hadley type cell in the meridional 
sea surface temperature was imposed, 
in the south. 


of the disturbance were set equal to 


and D) and Equations (2.22) through (2.27) 





Zu, y aF _ E 
y 0 7 aya Pp MIE 00 (2.30) 
and 
2 
BR, or ¡2 1 8 E a 
oe be ir, tom, a y 
y 
ie = | (aE _ 33E 
une ato 
Ss 2f Ap | a o | R (F + cp)! (2-31) 
For simplicity the following definitions were made: 
=- 8_ 
Ae > 
= ne pg 
a 
OS 2f Ap (2.32) 


hat a steady state condition existed 


in the model was made, i.e. — ) = 0, yielding 
rs. g | oE 2E 
0 = z15 w,| 0: E (2:33) 
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Since S vas mot equal to zero, 








1 
E 3E _ 5 OF 
Iw, | = 0 and ee as (2.308 
and Baustıon (2.31) can be written 
2 2f 
ð F O _ 
59 77S SH, i - (F + cp = 0 (2.233) 
dy 
or 
225 2f_ 2f 
S3 2 = 53H, 2 KE- Pa + SoH, R ci = 0. (2.36) 


Since F represents the mean temperature throughout 


the layer, the zonal mean thermal wind can be defined as 


T =- Z, 





q = a (27) 


Equation (2.36) was solved analytically, using the boundary 
condition that UN = 0 at y = 0, VW. 


The solution is 


1/2 y 
un se 1 sinh[ (8,55) (y-7) 
F= s. Sr os sy? ea Te ; 177 = ; 
where 4 5 2 
ar = 
= , S, ==-H, , 
oe c ay 4 5, 1 
2£, 
un: (2.38) 


This solution, along with the condition that A thru 
D equal zero and E equals 2F, was used as an initial state 
to expedite reaching a steady state. The initial mean 


thermal wind obtained from (2.38) was 


1/2 wW 
_ ou A sien U 
5 cosh[(5,S,) 7] 5 


where UL = 0.00433 degrees per kilometer and the change in 


the thermal wind with time was examined. 
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The complete equations were integrated in time 
starting from this initial state, After 40 days, values of 
Dr were within two percent of the initial value. The mean 
meridional circulation taken from the numerical solution is 
shown in Figure 3a. It consisted of a single cell with 
northward flow at the upper level and southward flow in the 
lower level. The vertical motion field consisted of upward 
motion in the southern part and downward motion in the north. 
The central portion displayed weak and variable vertical 
motion patterns and were interpreted as being meaningless. 
It follows from the equations, tħat the quantities A 
through D will remain zero. The method by which vertical 
motions and the’mean meridional motions were calculated is 
described in Appendix B. 

A similar experiment was made with the same linear 
sea surface temperature profile as before but with the 
initial conditions given by Equation (2.29a). This initial 
condition primarily differs from the previous one in that 
it contains a finite amplitude disturbance with wave number 
k in the barotropic part of the motion (A # 0). This in- 
tegration required more computer time to reach a quasi- 
steady state. A three cell structure as on in Figure 3b 
developed after 230 days. The northward and southward 
motions were stronger and the vertical motion patterns 


were more intense and well defined. 
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72 | 58 | 45 31 18 


Latitude 
cure 5.222, 7 Meridional circulation arrived at in Hadley 
aesemesexperiment. b., Meridional circulation of atmos- 


phere with disturbances included. 


A quantitative comparison was made of velocities 
obtained in the two experiments. Table 1 lists the values 
DER a. 
of the upper level meridional velocity in cm sec and the 


ya . -1 
mean zonal velocities in m sec at levels one and three. 
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MEAN MERIDIONAL VE- LOMAS eV IOC Lay 


LATITUDE JLOCITY UPPER LEVEL | ]EyEL 7 LEVEL 3 
HAD HAD DIST. 

72.0 0.0 0.0 O AO 0.0 0.0 
68.4 2.9 3.4 20.5 | 6.6 8.0 Cl 
64.8 3.5 1209 18.9 | 6.0 6.8 = 
61.2 3.6 16.0 oo 5.0 6.8 = 
5758 A oa | ee a, 5.1 
54.0 3.6 S53 MS 6.8 r, 
50.4 3.6 D7 19.3 |34.9 6.7 15.6 
46.8 3.6 ee 19.3 [36.8 6.8 16.6 
43.2 3.6 28.8 1 e ro 14.6 
39.6 3.6 DEE To 6.8 9.8 
36.0 3.5 -8.6 MOS 6.7 3.7 
32.4 3.6 6.4 Op? | |) pee 6.8 OS 
28.8 3.5 10.1 or e7 6.8 | -0.8 
25. > DS 3] 172099 1380 6.8 07 
21.6 2.9 1.4 OS 8.0 1.8 
18.0 o 0.0 OO 0.0 0.0 


Table 1. Comparison values of meridional and zonal velo- 
cities from Hadley experiment and disturbed atmosphere 
experiment. Positive values are northward or eastward_ 
flow. Upper level is 500-0 mb layer. Units are m sec 
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2 Energy Examination 


Equations for the components of potential and kine- 
tic energy and energy transformations in terms of the Four- 
ler coefficients are as follows: (The derivation of these 


equations is shown in Appendix C.) 


N 





pe dy=b s Fay, (2.40) 
2 
fe ay=# 5 (e ey) ake (2.41) 
_ 2 IE, 2 oF) 2 
re Cs.) 1 dy , (2.42) 
ae dA, 2 OB, 2 oC, 2 9D, 2 
dy > J I (ay) + oe T (Sy? ih (sy? 
+ k? (A? + B? + CÊ + D?)] dy, (2,43) 
Er = 
G = Hea A yr | dy, (2.44) 
-_r fk 42 £ oF 
> =-/[ [5 U (BC AD) T dy (2 245 
2 2 2 2 
c, 1 os EB _ A, (22D _ „ac, 
Y Jy? Jy? dy? dy? 
3 2 2 2 2 
- AZ De oe Ct ST l ay, (2.46) 
y dy? dy? dy? dy? 
ce pg Dre 2 
DT 5 4 (ok) EB) Le = ck") 
y 
2 
Ze prê) dy , (2.47) 
dy 
Ic). _ F lee (ae 
D 15 eC, [(4-c) (1,1 (A-2C) 


9 = ðC JA 
+ ay LIV, lA = DAD re) 


(15,12 a-20m + mie? - 21) flay, caer 
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°F 


= Lg meee 
De Ex] 4 (pk) (56) 9 LE m. hay, (2.49) 
u. 00 y) Lily SEM JE 
D = ri ec, [E-P 117,12 51] favs (2-50) 
ae 
coo -f [75 u Fl dy (ESAS 
and 
Zn 
Ca = -f E (wot aF wD) ] dy. (22525 


The relationship of these quantities can be expressed by 


the following equations: 


d 

— == o + 

ge PY IA A (25) 
a ec = 0 , eos) 
ot e p e j 
az, = C= pL, = D (2.55) 
ot e e k ei es ; 
and 

en TTK dy = -C +C, -D - D (2.56) 
ot Z Z k zi Zs ’ " 


T£ Equations (2.53) to (2.56) are summed, the following 


equation is obtained: 


0 = 
ae S (P+K) dy = Cc, Do; Do DE D (257) 


and shows that the time rate of change of total energy in 

the system is equal to that generated minus that dissipated. 
The quantity a represents the conversion of mean 

potential energy to eddy potential energy that is accon- 

plished by the horizontal transport of sensible heat by 

the eddies. The temperature in the layer one-three is pro- 


Portional to Ve and k(BC - AD) r is proportional to v T: 
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Qu The quantity k(BC - AD) represents v'T' which is the 


dy" 

transfer of heat across a latitude circle, and since F 
tends to be negative, k(BC - AD) must be positive to convert 
Eto P . 

Z e 


The Term Ea represents the conversion of eddy poten- 
tial energy to eddy kinetic energy by the vertical circula- 
tion in the zonal plane. The conversion depends on the 
correlation between w' and a' (a' is a departure from a 
mean specific volume of air). Since w'a' is proportional 


to On! > w, and C, and w, and D must be negatively corre- 


C D 


lated to convert e to Kae 


The quantity C, represents the conversion of eddy 


k 
kinetic energy to zonal kinetic energy and is dependent on 
the correlation between the mean zonal wind and the meri- 


dional convergence of the eddy momentum transport. The 


, =O 
guantity u ~ (u'v') is proportional to 


oy 
2 2 2 2 
ee ee ee obig 
yoy dy dy dy 
2 2 2 2 
+ kn (a2 - p24 - pot + cs E 
yay dy dy 9 y 
Since on and cla tend to be negative A) must be 
dy 9 y ” Oy 


positive in order to have conversion of E to K_. 

The term C, represents conversion of zonal poten- 
tial energy to zonal kinetic energy by the mean meridional 
circulation. It depends on the correlation between w and 
o, which is proportional to WO Pa This conversion is small 


nina tute ( Cort, 1964) and requires a positive correlation 


to convert K OS 
Z Z 
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Pesq ante Lty > is the generation of zonal potential 
energy through diabatic process. In this model the heating 
is a function of the air-sea temperature difference. The 
correlation of T and O, or F and Q, must be positive to 
generate zonal potential energy. This means that heating 
of relatively warm air and cooling of relatively cool air 
occurred in the model. 

The dissipation terms Doi? Dog? Den and De provide 
an energy loss mechanism for eddy and zonal kinetic energies. 
The second subscript, s or i, refers to dissipation by sur- 
face or internal friction respectively. 

An interesting way to examine the system energy is 
by considering an energy flow diagram. Figure 4 depicts 
swch a diagram. The boxes contain time rate of change of 
energy, and the arrows show the direction of the flow of 
Nery if thé transformation is positive. 

The model equations were integrated until a statis- 
tically steady state was reached. Energy calculations were 
made and instantaneous values for the quantities just dis- 
cussed are shown in Figure 5. 

The diabatic heating from the sea surface creates 
zonal potential energy that amounts to 1.27 watt N The 
conversion from zonal potential energy into its eddy 
counterpart occurs at the rate of 1.32 watt a The eddy 


potential energy is converted to eddy kinetic energy through 


the eddy circulation at a rate of 1.34 watt a The 
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Figure 4. Energy flow diagram. Flow is in the direction 
of the arrows if the transformation is positive. 
conversion of eddy kinetic to zonal kinetic energy, which 
maintains the zonal currents against frictional dissipation 
320,63 watt ag ihe mci ssipatihon af Teddy kinetic energy 
by surface friction is 0.79 watt na This is substantial, 
larger than dissipation by internal friction and the dissi- 
Parion of zonal kinetic energy by surface friction. 

It should be noted that the change of potential 
energy, both zonal and eddy, was very small while the change 


in kinetic energy was substantial. The atmospheric model 
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m 


605 7 


Figure 5. Flow diagram of the model atmosphere energy. 
Energy and energy transformation units are in watt m-2 
lg "ere cm’ sec’). 


had reached a statistically steady state at the time the 
energy computations were made, but the energy continued to 
fluctuate about a mean value. This undoubtedly influenced 
these calculations. 

An exact balance was not achieved throughout the 
entire system. It is believed that these very small dis- 
drepancies were due to truncation or roundoff errors and 
should not be considered a serious defect in the calcula- 


tions. Figure 6 is included to show that these values are 
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Figure 6. Flow diagram of atmospheric energy as con- 
structed by Oort(1964). Energy units are in 10° joule 


2 . e . — 
m Senergy transformation units are in watt m 


not only qualitatively correct but also quantitatively 
acceptable. Figure 6 shows the atmospheric energy flow 
diagram as constructed by Oort (1964). The values en- 
closed in the boxes are no longer time rates of change of 
energy but the amount of energy present. 

It is satisfying to note that the transformations 
in the model atmosphere compare favorably with those of 
Oort. The conversion os and C, agree very well while Co 


and Ce lie just outside the extreme values as proposed by 
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Oort. The dissipation of eddy kinetic energy is much small- 
ersthan Oort s. Since his value itself is considered too 
small, this points out a weakness in the dissipative mechan- 
ism of the model. This may be due to the lack of horizontal 
diffusion terms in the model equations. It was decided that 
diffusion was unnecessary in the model and this decision may 
be reflected in the energy dissipation. Also, energy loss 
through nonlinear interaction between waves has been elimi- 


nated in this model by the single wave formulation. 
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Dres ER OCEAN MODEL 


DES CREFPTION AND EQUATIONS OF THE MODEL 

The ocean model is similar to that studied by Bryan 
(1969) and es Choy Taye tt is ada baroclinic model based 
on the primitive equations in which steady state solutions 
are sought by direct numerical integrations of an initial 
value problem. The circulation produced in the ocean-only 
model is primarily driven by surface forces that simulate 
wind stresses, and horizontal pressure gradient forces due 
to internal stratification. The horizontal stress exerted 
on the ocean surface by large scale wind systems causes the 
surface water to move as an Ekman layer and drift to the 
right of the stress in the Northern Hemisphere. This causes 
the curl of the wind stress to become important as it causes 
convergence or divergence of the Ekman layer transport and 
corresponding vertical velocities below the Ekman layer. 

Heating in low Latitudes and cooling in high latitudes 
causes a thermal circulation to develop. A thermal current 
is set up where there are easterlies in the south and west- 
erlies in the north. These currents converging into the 
continental boundaries contribute to the resulting circula- 
tion, 

The rectangular ocean basin is everywhere 1700 meters 
deep and is 4800 km wide and 6000 km in length. To coin- 
cide with the atmospheric domain the basin extends from 18 


degrees to 72 degrees north (Figure 2). As in the model 
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atmosphere, the ß-plane approximation is used. Veronis 
(1963a,b) has shown that the B-plane approximation is satis- 
factory for the intermediate scale motions in middle lati- 
tudes. This formulation is also used to maintain consistency 
with the atmosphere and prepare for future joining of the 

two models. The -plane approximation in both models elimi- 
nates the need for mapping and reduces the possibility of 
confusion that might arise if an ocean model in spherical 
coordinates were joined with an atmospheric model using car- 


tesian geometry. 


2 
Fe A NM 500 
nn en sy eee ne ee: Wu. 900 
da 

Tuy 1300 
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Figure 7. Vertical structure of the ocean model. 
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The vertical structure of the three-level model, shown 
in Figure /, is changed considerably from that described by 
Haney. The three levels are placed at 100, 500, and 1300 
meters below an undisturbed surface. Temperature and hori- 
zontal velocities are predicted at these levels while the 
vertical velocity is diagnostically computed at the inter- 
mediate levels. The vertical layering is defined in such a 
way that there is a finer resolution near the surface than 
near the bottom. 

The Boussinesq approximation (variations of density are 
neglected except in terms involving gravitational accelera- 
tion) and the hydrostatic approximation are applied. In 
the equation of state, the dependence of density on salinity 
is neglected. Density is assumed to be a linear function 
of temperature only. 


me governing equations are 


du _ _ i op i 
7 5 a + yí + E, 9 (3.2) 
O O 
dv 1 op 1 
— = - — - uf + — F ; (332 
dt Po oy Pa y 
OB AN 
T Dg , (3030 
a OW 
V «Y + Ir = 0 , (34) 
Z 
aio 
dt oC ° 2 
O 
and 
p = p 11 - a(T - Tod: (3:60) 
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where F = eastward frictional force per unit volume, 
E zemeortebward frictional force per unit volume, 
and Q = heating rate per unit volume. 


The quantities Fo N and Q represent the assumed con- 
tributions of sub-grid scale motions to the exchange of 
momentum and heat in the ocean. The internal redistribu- 
Eons accomplished by utilizing a linear lateral diffusion 
oca and Momentum. Similarly, the vertical mixing of 
these two quantities by sub-grid scale eddies is specified 
by a linear vertical diffusion coefficient. The vertical 
mixing of temperature is enhanced by the use of an instan- 
taneous adjustment mechanism to eliminate unstable lapse 
rates, 

For the initial experiments with the ocean model a con- 
stant wind stress was imposed. Figure 8 shows the mean 
zonal stress that was used in this phase, The data was 
taken from Hellerman (1967). 

Following Haney (1971la) the downward heat flux at the 


ocean surface is expressed as 


Koz) HET, - T,,,)1/0,¢ « Gap 


The quantity H contains the net downward flux of solar 
radiation across the ocean surface, minus the upward flux 
of longwave radiation and latent heat from an ocean surface 
at a temperature egual to I,» the prescribed atmospheric 
equilibrium temperature. The second term contains the up- 


ward longwave radiation and sensible and latent heat. The 
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coefficient H, was found by Haney to vary with latitude by 


1 
about 20 percent. It was set equal to a constant in this 


work. 


72 


Latitude 





Stress 


Figure 8. Mean zonal stress used in the ocean model 


experiments. The data is from Hellerman MS Es 
dy nes cma . 


B. COMPUTATIONAL PROCEDURES 

Since the ocean bottom Hoer hate = Oise the boundary 
Condition on the vertical velocity there, Following the 
method used by Bryan and Cox (1967), Haney (1971la) and 


others, the ocean surface was considered a balanced surface 
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at which w = 0. The height of the surface was not obtained 
from the continuity equation so external gravity waves were 
eliminated, which permitted a much larger time step to be 
used. However, the technique also removed the vertical 
mean divergence from all scales of motion, including Rossby 
waves, and introduced an artificial constraint on the phase 
speed of these waves (Gates 1968). The necessity for the 
longer time step to accomplish long-term integrations made 
this technique mandatory. This filtering technique does 
not affect the thermally driven part of the motion. 

A new set of equations was formed by subtracting the 
vertically integrated form of Equations (3.1) and (3.2) from 
Equations (3.1) and (3.2) respectively. After eliminating 
the pressure deviation in favor of the temperature by the 
hydrostatic equation and the equation of state, equations 


for the time rate of change of the vertical shear current 


(u'v') were obtained. These equations have the form 
ay (3.8) 
ot za ' 
and 
9 ' 
7, (v ) = F,(u,v,T) , (3.9) 
where 
Del Rey ze toe. yu Zk) = ED (3.10) 
and 
ee E a E) = 9(x2,y,T), ES 
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and u and v are the vertical mean motions discussed shortly. 
Equations (3.8) and (3.9) were numerically solved to obtain 
new values of the vertical shear current. 

Oca y weduautonw fron which the vertical mean cur- 
rent was predicted, was formed by eliminating the pressure 
between Equations (3.1) and (3.2). The result was vertical- 


ly averaged and had the form 


O 
9 1 me 
SE iS Cdz = F,(u,v) (32125 


Are cC is the relative vorticity and H is the depth, of the 
basin. Vertical integration of the continuity equation 
(3.4) allowed the íntroduction of a stream function for the 
vertically integrated velocity since the vertical mean 
divergence was zero, The vertically averaged vorticity 
therefore had the form 

1 O 

if tdz = V*y (3.13) 


-H 
and Eqüation (3.12) can be written 


ð = 
aes “y= F (u,v) . (3.14) 


A relaxation technique was used to solve Equation (3.14) 
for the time rate of change of Y. The quantity Y was then 
set to zero on the boundaries, A time integration was then 
performed to obtain the new value of the stream function 


and then the vertically averaged velocity components from 


O 
“wi udz = - oY = u 


H dy 
-H 


and 1 O 
q f vdz 


zu 


I 


(3.2187) 


Ae 
IH 
<| 


52 





The total current was then obtained by adding the verti- 
cal shear components predicted by Equations (3.8) and (3.9) 
to the vertically averaged components predicted by (3.14). 

To obtain new values of the temperature, Equation (3.5) 
was written as 


m 
Ben 


(ee) (3,10) 
and integrated in time to obtain T. The quantities p and w 
yereithen diagnostically calculated from (3.6) and (3.4) 
\ 

respectively. 

The boundary conditions on momentum and heat at the 
walls were those of zero slip and perfect insulation. The 
ocean bottom was also considered to be perfectly insulated. 


As a simplification of the model, bottom friction was neg- 


lected and free slip was allowed. 


C. HEATING AND FRICTION 

It was assumed that all solar radiation penetrating the 
ocean's surface was absorbed in the top layer which was 
300 meters thick. Since the side walls and bottom of the 
basin were insulated, the only source of heat was at the 
surface and was calculated from the thermal boundary condi- 
tion (Equation 3.7). A parameterization scheme by which 
the internal redistribution of heat by sub-grid scale 
phenomena was represented by a simple linear eddy diffusion 
of heat and an instantaneous convective adjustment mecha- 
nism was used in the model. The right hand side of Equation 


Ocean be vritten 
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ER 2 8 dt AST) 6 , 
o C Ay VT t yz Ko) T A t c (32179 


O 
Bryan and Cox (1967) have shown that the lateral mixing 


coefficient, A for the apparent temperature should be 


m: 


less than that for momentum given by A This is due to 


M’ 
the larger scales of turbulent motion in the ocean having a 
tendency to take place along isentropic surfaces rather than 
demni Cetly horizontal surfaces. The lateral mixing of momen- 
tum and temperature in the model is somewhat dictated by 

mie resolucton of the numerical grid. The coefficient 
Should be sufficiently small so as not to obscure the 
lateral transfer of the quantities. A minimum value of 

both Ay and Ay is needed to satisfy the zero slip boundary 
conditions, and to simulate the horizontal eddy transport 

of heat and momentum by scales of motion too small to be 
resolved by the numerical grid. In this phase of the work, 
coefficients were held constant with 


8 2. al 


A 2S LO Tem Sec 


H 


and 


8 2 ai 


A 1:0 x 10 cm sec 


M 
Megverte cecal diffusion coefficient, K, differs in that 
the exact physical processes by which vertical exchange 
takes place within the thermocline is unknown. Following 
Bryan (1969) a value of 1.5 ae eee was used throughout 


this study.: 
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The convective adjustment mechanism (oT). is difficult 
to state mathematically. This device is a means by which 
the temperature structure in the model was forced to remain 
stable. After the convective and diffusive heating terms 
were added, the vertical temperature profile was examined 
for unstable lapse rates. If the temperature field resul- 
ted in an unstable lapse rate, the temperatures were in- 
stantaneously adjusted to the vertical mean temperature of 
those unstable layers. When this first adjustment was com- 
pleted, the temperature field was again examined to deter- 
mine if any other layers had been made unstable. The 
adjustment process was repeated until all layers were 
gravitationally stable., 

The quantity a was evaluated as follows: 


oz 


oT a 
az u IH, i H (T ec 4 





oT E - T 
Koz) k=15 = Ze, , 
=] 29 
= T = T 
oT 2 3 
Se) a al = a 
OZ k=25 2, Zo 
and 
aL p 


The last equality imposed the condition that the ocean 
floor was insulated since the vertical temperature gradient 
vanishes. 

The vertical eddy diffusion of heat at level one, for 


example, was then defined to be 
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K(T,- T,) 


as a De) 1/0 56 o (2,7 zo) 


3z (Koz? 1 1 


Ba Er (32,198, 


The heating term due to lateral eddy diffusion of heat in 


(3.16) may be written as 


2al ð ,oT ae oT 
A G t ay ay?) : (3.20) 
oT oT . * * 
The terms SE and dy were defined at points one-half grid 


distance to the east and north, respectively, of points at 


which T was stored (see Figure 9). The boundary conditions 


were 

nr = das > 

2 2 

and 

oT oT 

ee, + = == dx i-i - 5 (3.29 
where i = 5 is a point one-half grid distance west of the 
western boundary of the grid(located at i=1) and A one 2 
is a point one-half grid distance east of the eastern 
boundary (located at i=i ax)’ These conditions force ST 


to pass through zero at the eastern and western boundaries 
thereby producing thermal insulation. Similar boundary 
Bine oT 
condition on ye at the northern and southern boundaries 
ensure that the normal temperature gradient vanished at 
the walls, making them perfect insulators. 
E E 


Ihesfreıetıonal terms a and =a Ins Equatsons (3.1) Mand 


(3.2) are composed of two parts. These are the sub-grid 
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scale vertical eddy stress and the sub-grid scale lateral 





Mer uslon OL momentum. These can be expressed as 

F 1 OT 

= + A. Vu 
p Oz M 

O O 
and 
F 1 oT 
Po pP, o zZ M 


: oT ee. . 
As with the term KT» the quantities > and > were defined 


at the half integer vertical levels. The vertical eddy 


stress was prescribed at the surface and defined to be 


mero at the bottom. 


level as 


x, 1 = 

wu: 
"x, we 
Do 27 
es oe 
7 
es e 
Po = 


ana the meridional 


ty 
Dy k= 5 i 
ty 
o, k=15 i 
Usp 
0 k= 25 i 


tine Tonal stress vas defined at each 


prescribed surface stress, 


uS) 
E j 
(2, Z5 
(ua u) 
2 3 
0 : (3-23) 


Stress vas defined as 


Do. 

(v4 a 
Ferse) ® 
$ 2 

| (v, - V3) 
ar) ° 
2 3 
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= 0, (3.24) 


meus the prescribed stress was purely zonal. 
The x component of the term involving the lateral eddy 


diffusion of momentum can be expressed as 


ne CR U ce on 
The quantities <4 and S were defined one-half grid distance 


to the west and south respectively of the velocity storage 


points. The boundary conditions were 
as = ee 
dx’ i=l ex ” 
du 2u 
ii Ax i=1 neS ? 
max ax 2 
du 2u 
en 2 L eae 1 
ya Ly) j=15 ° 
and 
du 2u 
es Ze 0 es (3020) 
ey II max EN JA max 2 


These conditions impose the zero slip boundary condition 

on the velocity at the north and south walls and the no 
flux condition at the east and west walls, Similar defini- 
tions for the meridional components impose the zero slip 
condition on the east and west walls and the no flux condi- 


ron on the north and south walls. 


DESEEN TTESDIFEFERENCE TECHNIQUES 
A space differencing scheme that utilized a centered 
differencing with a staggered grid was used in the integra- 


tion of the primitive equations. When a centered space 
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scheme is used with an unstaggered grid, a computational 
mode in space arises that is analogous to the computational 
Bimesmode that occeurs with a centered time differencing 
(leap frog) scheme. The use of the staggered grid not only 
prevents the appearance of the computational mode in space 
but is responsible for considerable savings in computing 
time because variables are retained at alternate grid 
points. Arakawa (1966) has shown that maintaining integral 
constraints on certain quadratic quantities of physical im- 
portance is sufficient to guarantee nonlinear computational 
stability. The space differencing used in this model is a 
Simplified, energy conserving version of the scheme used in 
the UCLA general circulation model, as documented by Gates 
Geal (1971): 

Total energy is conserved in this model, The boundary 
conditions of zero normal flow through the side walls, and 
zero vertical velocity at the bottom, guaranteed that the 
following integral constraints were satisfied by the differ- 
ential equations. 


Kinetic Energy 


d al 2 2 
Ae | 2 Post * Vv dav = Jit-esw + p¿(uF, + vF )Jdv. (3.27) 
Vv v 
trotential Energy 
d- f og(+n)av = ff leew = eet i) av. (3.28) 
yv vV 


The potential energy was taken to be zero at the bottom. 


3:9 





ME added to (3.28), a conservation equation 


for the total energy in the basin is derived and is given 
by 


d 1 2 2 p 
76 p lu ty) e pg(z+H)dyv = 


V 
Jo, tor, + E? - g(z+H)®2]dv. (3.29) 


Beenlterditserence forms of the prediction Equations 
MS (3.9), (3.12) and (3.16) are shown in Appendix D. 
Haney (197la) has shown that these difference equations 
conserve total energy. 

Figure 9 shows the staggered grid network. The velo- 
city components, u and v, are located one-half grid incre- 
ment east and one-half grid increment north of the points 
where corresponding scalar variables T, p and UY are re- 
tained. The grid distance is 200 km in both the x and y 


G@irections. 





Figure 9. Location of variables in the horizontal 
Zeapgsered grid. 


60 





Kmerneseimerdirrereneinse scheme, various terms in the 
rar lonssoftswotrlomsewere treated differently. The pressure 
gradient term was evaluated using a centered difference. 
The Coriolis term was evaluated by a trapezoidal implicit 
scheme, Prd the friction term was evaluated using a forward 
time step to ensure linear computational stability. 

The trapezoidal implicit scheme for the Coriolis acce- 
iteration reduces to a form 


n+l n-1 n+l 


u - u = 2At[aif u + (l-a)if a 


}. 2720) 
It is easy to show that this scheme is linearly computation- 
ally stable for all At if $< a <1. In this model, 

Mee O55. 

To supress computational modes in time that can arise 
when a leap-frog scheme is used exclusively in an oscillat- 
ing system, a Matsuno (Euler-Backward) step was used 
periodically. This two-step scheme damps the computational 
mode and permits the long term integration. 

The leap-frog scheme is conditionally stable and the 
maximum time step was determined by the highest frequency 
governed by this scheme. External gravity waves have been 
removed from the model which leaves internal gravity waves 
as the highest frequency waves present. The phase speed of 
these waves is g'H' where g' is a modified gravity equal 


to „a Bine ssmodelesen. 3.60 cn sa 


Po 


meters and Ay = 200 km, and the linear computational sta- 


, H' = 1200 


Bility criterion is 


61 





A diay 3°. 
an 


This was the time step used in this phase of the model. 


At 


4A 


E. INITIAL CONDITIONS 

A state of rest was imposed on the ocean for the start 
of the first integrations. A stratified temperature struc- 
ture was selected so that a minimum amount of computing time 
would be required to bring the ocean model to thermal equi- 


librium. Figure 10 shows the vertical temperature structure. 


| | | | | | ` 
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Page lO. ed )einitial temperature field at level one; 
b) Same at level two; c) Same at level three. 
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The northernmost temperature at each level is 4 degrees 
Centigrade and the vertical structure is isothermal in that 
column. The temperature at level one increased from 4C at 
ene northern boundary to 26C at the southern wall. At level 
two the temperature increases from 4C to 8.5C and at level 
three it is a constant 4C throughout. The stability varies 
from isothermal in the northernmost column to very stable 
in the southernmost column. No zonal variations were ini- 


erally imposed. 
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IAE COMBINED AIR-SEA MODEL 


A combined model was developed by joining the atmosphere 
and ocean models. The exchange of momentum and heat across 
the air-sea interface is the mechanism upon which the cou- 
pling of the two models was based. An exchange of moisture 
was not possible since the atmospheric model was "dry", but 
a crude attempt at parameterizing latent heat transport in 
the atmosphere was incorporated and will be discussed in 
the section on heat exchange. 

The physical aspects of the combined models remained 
the same as they were in the individual models. The spatial 
grid in the ocean, and time increments for each model, re- 
mained the same. The time marching process allowed the 
time integration to remain synchronized and avoided a dis- 
torted response of the ocean circulation to atmospheric 
fluctuations. The atmospheric solution was evaluated at 
grid points in the x-direction in the atmospheric model. 
This was done to construct a complete stress field from the 
atmospheric Fourier coefficients in order to calculate the 
air-sea momentum exchange. The ocean surface temperature 
was transformed to wave number space to calculate the air- 
sea heat exchange. The overall techniques of integration 
used in each separate model remained unchanged in the conm- 


bined model. 
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BZ ME SYNCHRONIZATION 

Meore problems that arose in early air-sea models, 
such as that developed by Manabe (1969a,b) and Bryan (1969) 
gas the impracticality of using a marching process that 
would keep the ocean and atmosphere synchronized. Primitive 
equation atmospheric models require time steps on the order 
of minutes and cannot be realistically integrated in sequence 
with ocean models that have time steps on the order of hours, 
Because the computing time requirements are prohibitive. 
Techniques have been developed to simulate time synchroniza- 
tion but little success has been achieved in long-term one- 
to-one integrations. 

Since the present atmospheric model is quasi-geostrophic, 
a time step of one-half hour was utilized. The ocean time 
step was nine hours as discussed previously. The resulting 
ratio of atmospheric time steps to ocean time steps was 
erehteen to one, The marching process consisted of 18 
thirty-minute time steps in the atmosphere followed by one 
nine hour time step in the ocean. 

Conservation of momentum and heat in the exchange pro- 
eessSduüuring the integration was crucial to the development 
oMa realistic interaction model. Therefore, the quantities 
Such as the total surface stress and heating calculated dur- 
wne chen S atmospheric time steps had to equal the quantity 
as it was applied during the one ocean time step. The 
Dernodiby whach this was accomplished will be discussed in 


ehe mnmext two sections. 
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B. MOMENTUM EXCHANGE 

The surface stress discussed earlier and computed from 
Equation (2.3) was used as an upper boundary condition for 
the ocean. The calculation of this stress was based on the 
assumption that there was no motion in the ocean. That 
assumption is also used in the stress formulation in the 
combined model. Since the velocity of ocean currents is 
much Smaller than the velocity of air flow, the effect of 
ocean currents on momentum exchange can be neglected. 
Manabe (1969b) used the same argument in justifying his use 
of a similar formulation of the stress. Simplicity and 
computing requirements also entered into the decision to 
use this formulation here. 

Starting from a prescribed initial state, the atmos- 
pheric equations were integrated over 18 time steps. The 
Surface stress was calculated each time and an average was 
computed to be used for the ocean integration. A general 
numerical procedure for carrying out the computation of 


the time mean surface stress is given by 


1 NM . 
er = ny y= CT) ijn] CELS 
and 1 NM 
Sir, = nl“ ne (4.2) 
EN Ce N correo sphere 


andi and j are horizontal grid indices having their usual 


directions, 


66 





Pee wtneomsiuErace Stress Components in the atmospheric 
model, e and Be: were not calculated directly, terms had 
to be collected and values for the stresses recovered at 
each grid point. The stresses were computed in three equa- 
mos. (2.22), (2.23) and (2.26). 


The following definitions were made: 


lii 


TOT: Ty), Conter buc ions of the A disturbance to the 


zonal and meridional components of stress 


computed in the A equation. 


T In Ty) y Coni buc ons Of the Bo disturbance to the 


zonal and meridional components of stress 


computed in the B equation. 


N} 


Ty O Concributions of the mean field to the 
x’E NOTE 

zonal and meridional components of stress 
computed in the E equation. 

To reconstruct zonal and meridional stress fields from 


these quantities, the terms must be multiplied by trigono- 


metele Tunetions as follows: 


aS i = Er + zen cos kx, + TL)» sin kx; (4:3) 
and 
(tay = Un: + GR cos kx, + e sin kx, . (4.4) 


The locations at which the atmospheric surface stresses 
Mone nmeateuilared coincide with the points in the ocean at 
which the temperature and stream function were calculated. 
mmecenuladnd yo were Calculated at the staggered or half-grid 


points in the ocean system (see Figure 11), it was necessary 
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to average the stresses calculated in (4.3) and (4.4) before 


apply ing the stress to the ocean. The result is 
Be area). 
X 1j 4 XSi] KLS Mee jc oe Peco 
ES (4.5) 
oy 3 ty A aaa) 
(4.6) 


SEXXEanNdgSSTXy were then calculated from (T )" and Ca yoy 
x ij ea 


for application to the ocean. This procedure of averaging 
over the number of atmospheric integrations, and applying 
the stress at the proper position was continued throughout 


Emer duration of the long term integrations. 


C. HEAT EXCHANGE 

To satisfactorily conserve heat in the air-sea exchange 
process, a time averaging procedure similar to that used in 
the momentum exchange was utilized. The equations for the 
heat flux in both atmosphere and ocean must guarantee that 
the total heat lost by the ocean equals that gained by the 
atmosphere. 

The atmospheric heating was previously discussed and 
eeen Eguation (221): In the coupled model, that for- 
mulation was expanded to include the net radiative heating 
of the atmosphere and a crude specification of the heat 
meteased through precipitation processes. The heat flux 
into the atmosphere was given by 


E eet HCE = 


ea 4? Wook, 
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Atmosphere Level 4 





_e Ocean Level 1 





<| 7 
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Figure 11. Grid alignment of atmosphere and ocean to show 
Bon coincidence of stress points in atmosphere and velo- 
city storage points in the ocean. 
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"Here Ro= net flux of radiation into the atmosphere. 
acdsee = net flux of heat into the atmosphere due to 
wee a pitation. The specification of R and L * P are des- 
cribed below. 

The ocean heat flux was also modified from the form 
used in the ocean only model. The necessity for modifica- 
tion arose from the fact that in the ocean-only model (like 
the atmosphere-only model) there was no requirement for an 
exact accounting of heat exchange. The ocean could give 
Hpeany amount of heat into the non-existent atmosphere with 
no adverse effects. Similarly, when heat was required, 
there was an endless supply available from the same non- 
existent atmosphere. In the combined model, the loss of 
heat from the ocean by evaporation was a specified function 
ey only. The quantity H_ in Equation (3.7) was divided 
into a term representing the net downward flux of short 


and long wave radiation, S minus the upward flux of 


N?’ 
Patent neat, L-E. The equation for downward oceanic heat 


Flux, 05 was therefore written as 


Oe ae See eh H, CES u (428) 


O N ea 4? 


where Sy znel dovnvard radiative heat flux 


and L'E net downward latent heat flux. 


Figure 12 gives a schematic representation of the heat flux 
into the combined model and the calculated exchange 

ET ea Ta) nemar rows Originating or terminating out- 
Side the model atmosphere or ocean indicate externally 


specified heat fluxes. 
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Ocean 
LE 
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meepure 12. Schematic representation of simplified heat 
flux into combined model. See text for definition of 
symbols. 

Annual mean values were used in the radiation terms in 
Pera Equations (4.7) and (4.8). Wetherald and Manabe 
(1972) recently reported on the effects of seasonal varia- 
tion of the solar radiation on the joint ocean-atmosphere 
model developed by Manabe and Bryan. Their model atmos- 
phere showed a marked warming in high latitude. Although 
a similar result might be expected in this model, varying 
the solar radiation should be the subject of a separate 


investigation and will not be included in this work. 
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From Equations (4.7) and (4.8), the total heat flux into 
the ocean-atmosphere system, Q, dr 05 LS given by 


Ono RL E Lp. (4.9) 


When Equation (4.9) was integrated over the horizontal do- 
main of the model, the total prescribed heat flux into the 
combined system was obtained. Time integration of the 
model cannot approach a steady state unless this total pre- 
scribed heat flux into the system is zero. Also, the atmos- 
pheric circulation is insensitive to the value of the total 
heat flux because only an input of available potential 
energy impulses the circulation. Since each of the heating 
components on the right hand side of (4.9) represent essen- 
Ata Ty different physical effects, each term was prescribed 
to give a zero net heat flux when integrated over the hori- 
zontal domain. The meridional mean was removed from R and 
Sn and the quantities L*P and L'E were made sinusoidal 
about zero in the y direction. Observed data shows that 
the ocean receives radiation and the atmosphere loses it 

at all latitudes. Removing the mean and using the new 
quantities R' and SN simulates the role that the media 
play as both a source and sink of radiative heating. Re- 
moving the mean from both L*P and L*E was necessary to 
insure that the net prescribed heating in both the atmos- 
phere (given by (4.7)) and ocean (given by (4.8)) was 

zero. Thus, there is no net prescribed heat exchange be- 
tween the two media - all of the heat exchange is calcu- 


PiEcdminternaliyetrom the last term in (4.7) and (4.8) and 
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i Vemecmonemt Men ath-sed temperature difference. Denoting 
new quantities that have a zero meridional average by 


primes, Equations (4.7) and (4.8) become 


R' + LeP'+ H. (T = fp 


a i "sea n? AE 


and 


Q =5 '-L'E'-H CT Sa 


O n if sea aoe (4.11) 


These terms will be discussed in the next sections. 
Peaeeinecspnertc Radiataon, R', and Solar Radiation, Sn 

The observed annual mean net atmospheric radiative 
cooling plotted as a function of latitude is shown by the 
dashed line in Figure 13. The data is from London (1957). 
It shows a meridional average loss of heat due to radiative 
effects of 184 langleys A The dot-and-dashed line in 
this figure portrays the deviations from the mean radiation. 
As can be readily seen, the deviations from the mean are 
very small and are in all cases less than 10 langleys aaa 
These deviations were considered too small to be signifi- 
emand vere neglected. Therefore, R* = O for future 
Salculations, 

Figure 13 also shows the observed annual mean solar 
radiation minus the mean upward longwave radiation at the 
earth's surface, Sy: This quantity with its mean removed, 
Sn is shown by the dashed line in Figure 14. The values 
of this quantity used in the model are shown by the solid 
Moe o igure. A profile that was both linear and 


Symmetric was chosen and held constant in time. 


ES 





day | 


Langleys 


200 


100 





Oana in 
-100 
-200 
72 45 18 
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Higure l3. Net mean atmospheric radiation, R, (dashed 


lernte). devsjation (rom mean atmospheric radiation, R”, 
(dot-and-dashed line), annual mean net radiation arriving 
e ace nS a (solid line). Data from London (1957). 
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72 45 18 

Latitude 
Deu e Deviation from mean net radiation, S,,', (dashed 


line) as observed by London (1957); model annual mean net 
radiation, I: (solid line). 
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e a enea enat Terms, L°’E' and L+P! 

The parameterization of the latent heat processes 
could be accomplished in one of three ways. Since the at- 
mospheric model contained no moisture, a direct exchange 
of this quantity and the associated releases of heat could 
not be considered. One possible formulation was based on 
the assumption that the air directly above the sea surface 
Bee sarıratred so no evaporation existed (L*E'= 0). This 
formulation coupled with the assumption that no precipita- 
Bon existed thereby eliminated the terms L°E' and L*°P' 
from Equations (4.10) and (4.11). An Zope mmen al shoni 
term integration of the model equations was completed using 
this formulation. Results were unsatisfactory in that it 
became evident in this short test that an active air-sea 
model could not be developed with this scheme. The merid- 
ional gradient E the sea surface temperature became weak 
and the atmosphere consequently became baroclinically 
Saper Atter a short time this formulation was discarded. 

Aaformulation of L-E! and LP that was similar to 
the observed profiles (Figure 15) was used in a second 
short term experiment. The result was similar to the case 
Just described, only the meridional temperature gradient 
was weakened at a more rapid rate. The solar radiation 
absorbed at the ocean surface, Sy was insufficient to main- 


tain an adequate meridional sea surface temperature gradi- 


ent and the addition of latent heating in the north and 
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cooling in the south further reduced the atmosphere's 
meridional temperature gradient. The results of this test 
were also unsatisfactory. One suggestion for the reason 
that the ocean did not maintain a strong temperature gradi- 
ent even though the applied solar radiation was close to 
that observed in the real atmosphere is that the model can- 
not simulate the sharp separation of the western boundary 
enrrent from the coast and the resultant tight gradient of 
temperature in the mid-latitude regions. This was probably 
due to the degree of linearity of the solution in this area 
(Ser Section VI A 1) and resulted in an oceanic and atmos- 
pheric meridional temperature gradient that was too weak 
for baroclinic instability. The system as a whole needs 
the latent heat released in the tropics of the real atmos- 
phere. 

The results of the above two experiments indicated 
that the ocean-atmosphere model required a somewhat differ- 
ent partitioning of the heating than is typically observed 
in the real ocean-atmosphere system. The heating formula- 
tion adopted for the model was equivalent to assuming that 
evaporation occurred in the ocean part of the model and 
EateenRezmosrstüre immediately precipitated out at the same 
latitude in the atmosphere (L-E' = LeP'). This formulation 
is somewhat unrealistic since, over the latitudinal domain 
ne mode Mathe observed evaporation exceeds precipita- 
tion in lower latitudes while the reverse is true in the 


mad kenlatituüdes. Figure 15 shows curves of the latent 
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Heat flox due to evaporation (solid line), and precipita- 
wono dashed line), from Budyko (1962). The meridional 
average values have been removed from these values. The 
As oda both L'E! and L'P' that was used to speci- 
fy these fluxes in the model is shown by the dot-and-dashed 
meen in Figure 15. While this formulation is quantitatively 
in agreement with the evaporation curve, it is a poor ap- 
proximation to the precipitation related heat release. How- 
ever, the greatest precipitation is in the equatorial region 
lying south of the wall at 18N. If the hemispheric preci- 
pitation were mapped into the domain of the model, the pre- 
scribed latent heat flux due to precipitation would be more 
in agreement with observed values. This is considered to 
be a weakness in the heating formulation and is an area 
that could be improved on in future work with this model. 
Nevertheless, it was believed that once the heating gradi- 
ent in the atmosphere was sufficient to maintain it in a 
baroclinically unstable state, the resulting air-sea inter- 
actions would be realistic. That is, the interactions 
depend on the characteristic atmospheric motions, independ- 
mie oi fthe external forcing that maintain those motions. 
Figure 16 shows the net prescribed atmospheric heat 
A dashed line) and the net prescribed ocean 
keat flux, Sy iC midi we dor to its final 
acceptance this form was tested in a short term experiment 
Nodos atistactorily contribute to the conditions 


teste foal wactive air-sea model, This fact coupled 
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writer iies (dot=and-dashed line). All have the meridional 
averages removed. 
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with the requirement for simplicity led to its use in the 
remainder of this work. 


cient ol el Reat Transfer, H} (T - T 
ee et 1s e a ——— 


a 


The quantity H, CI T,) will be referred to as 


sea 
Eher sensible heat. In the derivation of Hy by Haney i019 72s. 
it was shown that this quantity actually contained the net 
upward flux of longwave radiation and sensible and latent 
heat which occurs due to the air-sea temperature difference. 


Simce the parameter H. varied little, a constant value was . 


| 


assumed. As in the ocean-only case a value of 55 langleys 
day was chosen. 

As shown above, all of the heat exchange between 
the atmosphere and the ocean must be accomplished through 
the sensible heat term. It is in this term that a strict 
accounting of heat exchange had to be done. 


The sensible heat flux applied during an ocean time 


step was the average of the heat flux calculated during the 


18 previous atmospheric time steps. The equation utilized 
was 
1 NM 
ai — NM ai (ren o O ? Can 
where NM = (At) ocean/(At) atmos = 18. Because of the 


x-wave formulation used in the atmospheric model the quan- 
ELEY Mosa - T,) had to be calculated in several Steps. 
eue imne duat ions (2.17), (2.19) and (2.20) an equation 


for T; could De written: 


a = a SU cos kx, eee sin kx, ; (4.13) 


ol 





where 











2f oT 

(TF), = = E, + 250.3 - Soe Sz), 
Oat 

ee ic 

4 J R j 
and 

2f_ 

Sr = Le $ (4.14) 


To be consistent with the wave number representa- 
tion in the atmosphere and to facilitate coupling the 
models through a single wave number, a Fourier analysis was 
performed on the sea surface temperature. At each latitude, 
the sea surface temperature was analyzed in x-space to 
yield a mean zonal sea surface temperature and the Fourier 
coefficients of its first mode. This can be symbolized as 


ea i = SE Dee u a kx,.(4.15) 


The various modes of the sensible heat flux in the atmos- 
phere were calculated from 
BT oa y T,) E H ECTS, 


Pas ee 


+ (Teat - DR cos Kx. 


+ (T ea? - DE sin kx, ] (4.16) 
Muguvenremapo lated ın Equations (2.24), (2.25) and (2.27) 
during each time step. The right hand side of Equation 
(4.16) was averaged over the 18 atmospheric time steps for 
In corporaelonsın Equation (4.22). This in turn was applied 
imsthe oceanic heating Formulation. 


The model equations were integrated for a short 


time period (300 days) and zonally averaged values of the 
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Time momiecatewere calculated. The results are shown in 
Figure 1/ by the solid line. The dashed line shows the ob- 
served zonal mean sensible heat flux as presented by Budyko 
(1962). These were added to the imposed values of L-P' to 
obtain the zonally averaged net atmospheric heating. Values 
che Observed sensible heat flux estimated by Budyko were 
added to the observed values of the radiative and latent 
heat contributions, and after their meridional averages 
were removed, these new quantities were plotted in Figure 
18. The solid line shows the observed atmospheric heating 
over the ocean and the dashed line shows observed atmospher- 
ic heating over ocean and ae combined. The dot-and- 
dashed line depicts corresponding values of L*P' + Hı (T ea 
- T,) generated in the model. In the northernmost segment 
of the domain, the model closely resembles the observed 
values over the ocean. From a region near 45N to the 
southern boundary, the model is in closer agreement eh 
Ehe data from a land and ocean environment. 

There are several possible explanations of the 
Baer in the model atmosphere heat flux vor SON- The 
Horm kation Of the net solar radiation as shown in Figure 
13 favors a sea surface temperature structure that would 
be colder than the atmosphere in the higher latitudes and 
Warmer than the atmosphere in lower latitudes. This would 
lfeddsboenegativye contributions from the sensible heat term 


ate her latitudes. The curve dips sharply in the area 
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Mestarlux.zauerto sensible heat exchange. 


Figure 17. 
Model (solid line): 


Hositive is flux from the ocean; 
Observed by Budyko (dashed line). 
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Figure 18. Observed zonally averaged atmospheric heat flux 


over ocean (solid line); observed zonally averaged atmos- 
Pherie heat flux over land and ocean, combined (dashed line 
model zonally averaged atmospheric heat flux (dot-and-dashed 
Mie observed data from Budyko, 1962). 
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where the solar radiation and sensible heat passes from 
positive to negative values. 

The formulation of the latent heat released by pre- 
Papiro ia the atmospheric model contributes to héating 
meine south and cooling in the north. The sharp dip in 
the total atmospheric heating occurs at the latitude where 
the prescribed atmospheric heating becomes negative. 

The sensible heat term may contribute to this 
abrupt decrease in another way. Large positive values of 
this quantity are seen in some regions and large negative 
values in others. The large negative values are physically 
unreal because the real atmosphere does not give large 
amounts ef heat to the ocean through this term over large 
areas. These values are considered departures from a mean 
and can be interpreted as small heating where the values 
are negative and larger heating where they are positive. 
The prescribed heating of the combined model has already 
been shown to meridionally integrate to zero. There can 
be no net heat flux by the sensible heat term if a steady 
state condition is to be reached. If a different formulas 
tion for the prescribed quantities was chosen so that the 
net specified heat flux into the system was not zero, then 


the net heat flux by H, (T - T,) would have to be adjus- 


1 sea 
ted. In the present formulation, however, there is no need 


to change the mean temperature of the atmosphere and the 


interpretation just discussed is justified. 
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teure l shows the result of subtracting the sen- 
Sible heat term from prescribed radiaťion and latent heat 


contributions in the ocean model. ire cunve has the 


' 


general appearance of the Sy 


Curve in Figure 14.> But a 

amp appears in the same area as the "dip" in the atmos- 

phere curve and should be expected for the same reasons. 
The oceanic heating formulation that was arrived 


at empirically and through experimentation was utilized in 


the remainder of the numerical study. 
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Figure 19. Average of heating in the ocean model. 
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VeBEROCEDURES 


me FHASE I 
1. Variations of the Combined Model 

It became evident during the short term experiments 
discussed in the last section that the computing time re- 
quired to achieve a state of quasi-equilibrium was excessive. 
Bryan (1969) has shown through scale analysis that the ver- 
tical advective time scale is on the order of decades while 
the vertical diffusive time scale is on the order of cen- 
turies. Both of these time scales must be relevant in the 
integration of the coupled ocean and atmosphere over long 
periods of time. The relatively shallow ocean depth and 
stratified temperature structure specified as the initial 
conditions enhanced the approach to an equilibrium state. 

The combined model with a 200 km grid distance and 
time increments as previously described required approxi- 
mately one hour of computing time for every 100 days of 
Simulation on the IBM 360/67 computer. This would require 
Doo hours Of operating time for 100 years of simulation in 
this phase of the experiment, and this was prohibitive. 

Several ways to reduce computing time requirements 
were considered. These considerations involved the ocean 
model since the atmospheric model was formulated simply 
and could be integrated rapidly. In a combined model such 


as this, almost any change in either time increment or grid 
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distance in one part of the model would require an adjust- 
Peso zthatsauantity in the other portion of.the model, 

A time reduction could be realized by linearizing 
the horizontal equations of motion, as was done by Haney in 
the first phase of his ocean experiments. This change was 
incorporated into the model, along with the assumption that 
the vertical mean motion was steady so there was no require- 


ment to relax for Y the time rate of change of the stream 


t? 
function. It was this feature that produced the time reduc- 
Elon. This Simplified form of the combined model was 
tested in a several hundred day simulation and the model 
advanced 140 days after one hour of computing time. While 
this was a fair improvement in computing requirements it 
was not enough to warrant changing the model in this way. 
Changing the distance between grid points was con- 
sidered. The model previously described had a 200 km grid 
mestamece in both x and y requiring horizontal fields of 775 
points (25 by 31). If the grid distance was changed to 
moo Km in the y direction and remained 200 km in the x 
Mmianection, the computing time would be reduced. The reduc- 
tion was due to the fewer number of calculations to be made 
(now 400 grid points) and the doubling of the time incre- 
ment in the atmospheric portion of the model. For predic- 
Peon purposes the atmosphere had grid points only in the y 
mecena e this increment was doubled, the time step 
ceou lda lso be doubled and the solution still remain stable 


tomal innear sense. The time step in the ocean, however, 
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must remain the same since the grid distance in x was not 
changed. This version of the model was tested and produced 
a 320 day forecast in one hour. This was a substantial 
time savings but it was believed that an even more acceler- 
ated pace could be obtained. 

rne final modification to the original combined 
model was made by changing the grid distance in both the x 
and y directions to 400 km. When the time steps were 
doubled to one hour in the atmosphere and 18 hours in the 
ocean, computing time was reduced by a factor of about 
seven. A test showed that this version of the model pro- 
amecedean ooo day Simulation in one hour of computing time. 
This was a satisfactory improvement and this form of the 
model was utilized in the early stages of the experiment. 

eo arse Crid (400 km) Model 

Munk (1950) has shown that the development of a 
Vector (frictional) boundary current of several hundred 
kilometers is to be expected. It was assumed that the 
western boundary current was the width of the new grid 
distance, 400 km. It has been shown by Munk that the width 
che boundary current (L is proportional to the cube 
tooto the ralio of the lateral momentum diffusion coeffi- 
cient to the gradient of Coriolis parameter. This can be 


expressed as 
A 
MAS 
La {£ 7) : Corals) 


IT 





Mio ee western boundary width to the grid size, 


31/3 


Ll8 Z, is proportional to [A,/BCLx) In orderio 
keep this ratio the same in both the 400 km and 200 km grid 
models, and thereby have similar horizontal resolution in 


both models, the following must hold, 


A = 8 A : (OFZ) 
(400) M200) 
; 8 2 ll 
nee a value or 1x 10 cm sec was used for Ay with the 
A00 km grid, Au 0% Ho was used as an initial value in 
this model. Similarly A, was increased so that a value of 


H 

Ay = l x 10° nr Fels was used in the coarse grid model. 

The initial conditions discussed previously were 
modified to conform to the new 13 by 16 point grid and the 
long term integration was started. Because of the large 
horizontal eddy diffusion of heat, the meridional sea sur- 
face temperature gradient dropped from its initial value of 
22 degrees over the 6000 km length to 18 degrees in the 
first 400 days. To have an active air-sea interaction 
model the meridional temperature gradient in the ocean must 
be large enough to force strong meridional temperature 
gradients in the atmosphere so that a thermal wind larger 
arms Barsrequiredstorsbaroelinie instability can exist, 
Such a large temperature gradient could be maintained in 
the ocean only by reducing the diffusion coefficient. Hows 
w a haee oeil icient cannot be arbitrarily reduced because 


in the reduction, the western boundary width as given by 


(5.1) may become smaller than the grid size. 


92 





Meeps e reduction in the lateral diffusion co- 
efficient was made. The solution was examined for computa- 
tional stability after each step. At the end of the first 
400 days, the lateral diffusion coefficients were halved, 
the decrease in meridional temperature gradient ceased and 
the solution remained stable. At the end of 1200 days the 
eoefficients were again reduced and at the end of 7200 days 
they were reduced a final time to values of Ay = 1 x Oh 
and Ay id Xk 10° with no effect on the computational 
stability. It was evident, however, that with the 400 km 
grid, the western boundary current was only marginally re- 
solved. 

The model equations were integrated over the coarse 
wid ror a period of 20,200 days (~55 years). The volume 
average of temperature in the ocean basin was one parameter 
used to measure the approach to equilibrium. This water 
temperature had slowly cooled from an initial value of 
6.3076 degrees to 5.2017 degrees at day 20,200. In the 
early stage of adjustment the cooling was at a very rapid 
tate, 5.2 degrees per century, but by the end of this part 
of the experiment, the net cooling rate had reached zero. 

Another parameter monitored was the meridional 
gradient of sea surface temperature. This quantity had 
Stabilized at 21.6 degrees per 6000 km after day 15,000 


amd remalned so until the end of this time period. During 


this time the atmosphere remained very active. A thermal 
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wind, (U, - U,/2), with maximum values greater than 11m 
eae, which is well above the critical value for baroclinic 
ws tability, was present in the mid-portion of the domain. 
A strong upper level jet and weaker lower level jet were 
present at the same latitude. The state of the atmosphere 
and ocean was very encouraging and the data for day 20,200 
eaS used to initialize the fine grid (200 km) model to con- 
tinue the experiment. More discussion of the state of the 
ätmosphere-ocean System at that time is included in the 
section on discussion of results. 

3. Data Retention 

Poo con tinuine wleb the description of the pros 
cedure used with the fine grid model, a brief discussion 
cata retention must be included. In a project requiring 
ler se amounts of computing time, it is mandatory that a 
great deal of consideration be given to which data should 
De saved. One must decide not only what parameters to 
retain, but how often and in what form. It is virtually 
impossible to repeat an experiment to extract parameters 
that were not saved. 

The data obviously must be stored for analysis at 
the end of the experiment. Data must also be retained for 
esen restarting the model. It should not be expected 
Eratcsa Lone—-term integration of this nature could be com- 
cano e Computing time period. The equations must 
Bemivterratred and data stored periodically so that after a 


reasonable operating time, such as one hour, the model can 
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PemcotOppe@mandsrestarted at a later time using the last 
data in storage as initial values. 

Throughout the experiment with the coarse grid 
model, data were stored at 100 day intervals. The ocean 
variables saved were the stream function from which the 
vertical mean current was obtained, the zonal vertical 


Ehe amseurrene, u , the meridional vertical shear current, 


and the temperature, T, at all three levels. The at- 
mospheric variables saved were the Fourier coefficients, 
ome Hoe ihe retention of Such a large quantity of 
data posed a problem, but it was remedied through the judi- 
cious use of magnetic tapes. As one clan of the experi- 
ment was completed, the data that would not be needed until 
the final analysis was moved from disk pack temporary 
storage to a set of magnetic tapes where it was permanently 
stored. The necessity for this became even more evident in 
the fine grid model integrations when the fields were in- 
creased from 208 data points per level to 775 points per 
level. Saving this much data would have been impossible 
without a large storage capability. 

The last problem concerning data storage arose dur- 
mg the final part of the long term integration. After a 
Saree OreGguwasct-equilibrium was approached, five points were 
chosen for retention of daily data for the last six years 
of the experiment and were carefully examined by spectral 


techniques. Data from an ocean-only and an atmosphere-only 


experiment were also saved over the last six years. 
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ne Grid (200 km) Model 

The data for day 20,200 from the coarse grid model 
Near iso terpolated gn both x and y to obtain initial 
melds for the fine grid model. The lateral diffusion co- 
efficients for both momentum and heat were held at the 
final values as used in the coarse grid model. While these 
coefficients might have been reduced over a time period as 
they were in the coarse grid, this was not done in order to 
avoid the readjustment that would be required. A further 
enange in the diffusion coefficient would increase the ad- 
justment period. It has been shown by Manabe (1969a), 
Gates (1968) and others that the values chosen previously 
are appropriate for the grid size used, 

The equations were integrated over the fine grid 
model for another 6500 days with data retained every 100 
days. At that time, 73 years, the ocean-atmosphere model 
Was approaching an equilibrium state. The volume average 
of the temperature initially started to slowly decrease 
but after a few hundred days it commenced an increasing 
trend that it maintained throughout the remainder of the 
experiment. The warming trend began at the rate of 0.8 
degrees per century but this decreased to about 0.5 degrees 
Per century at /3 years. Compared to the warming rate ot 
dA degrees per century arrived at by the ocean in Manabe 
and Bryan's combined model, this was considered an accept- 
anent ime to Start the last six years of the combined 


model experiment. 
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pera at the 7/3rd year, a final 2200 day (6 
Metas) "Simulation was produced with the fine grid model. 
Tt was during this time that not only were all the data 
moconrdedeevery 100 days but for five selected points, the 
data were recorded every 1.125 days (every three oceanic 
time steps). These points were selected as being represen- 
tative of features that existed in various parts of the 


cean. Figure 20 shows the location of the five points, 


72°N 





6000 
- 18°N 
=: 4800 ——— 
Ergure 20, Location of points where selected data were 


Meco ded every 1.125 days. 
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E through E. Point A is the northernmost point in 
mie eastern “Sector where sea surface temperatures were 
warmer than the initial values and where strong sinking 
Roetone mad occurred. Point B is in the center of the domain 
where the cyclonic scale atmospheric disturbances had their 
mecatest eifect. Point C is in the western boundary current 
amd, if the model domain is likened to the North Atlantic 
Ocean, this point is close enough to an observation site 
used by Woods Hole Oceanographic Institution so that some 
comparisons can be made. Point D is an area in the south 
central portion of the domain and would be equivalent to 
warea near Bermuda. Point E is in the warm core of the 
ocean close to the southwestern corner of the domain. The 
data at these points were examined spectrally and the dis- 
cussion of this technique is included in a later section. 

At the conclusion of this part of the experiment, 
the combined model had been used for a 79 year integration 


of the equations, This completed phase I of this project. 


bere PHASE II 

This phase of the experiment consisted of integrating 
the equations of the ocean-only model for a six-year period 
uSing the fine grid. To simulate the coincidence of this 
emee od with the final six years of the combined model, 
data from the combined model were used for initialization. 
Initial fields of temperature and velocities were taken 


Boom day 26,700. 
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A fixed surface stress and an atmospheric temperature 
fadebo be imposed. Rather than use Climatology, the input 
ess and temperature were obtained from the combined 
model. In preparation for this phase of the experiment, 
the surface stress, and atmospheric surface air temperature 
were averaged over a 100 day period centered at day 26,700. 
These average values were then held fixed during the entire 
eme integration of the ocean-only model. 

Data were retained in exactly the same manner as in the 
combined model. Full fields of the ocean variables were 
woede very 100 days and data for the five special points 
were stored every 1.125 days. The variables were examined 


and comparisons were made with the results of Phase I. 


These results are discussed in a later part of this work. 


ass PHASE III 

The equations for the atmosphere-only model were inte- 
grated over the same six year time span with the fine grid. 
Initial values of A through F were taken from day 26,700 
of the combined model data and the 100 day average sea sur- 
face temperatures were imposed as a lower boundary condi- 
tion. Data were saved at the same times as in the first 
two phases. 

Galicia tiones tol energy and energy transformations were 
mader and are described later. The data from the five points 
denena azed spectrally and this too is presented in a 
te a cection. With the completion of this phase of the 


experiment, the long-term integrations were concluded. 
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DeSSEREVTEW OF COMPUTING TIME REQUIREMENTS 

Pmccewsehe  inteprationszwere condueted with a warietyzof 
Zumerırecar models, it is germain at this point to review the 
computing time requirements for each model. This review is 
Presented in tabular form in Table 2. 


LENGIH OF SIMULATION 
MODEL DESCRIPTION FOR 1 BR COMPUTING TIME 
























































Combined Ax=Ay=200km 101 days 
(Pine Grid) At =30 min 
atm 
At = 9 hours 
ocn 
Combined Ax=Ay=400km 666 days 
(Coarse Grid) At = ] hour 
atm 
At =18 hours 
ocn 
Ocean Only Ax=Ay=200km 308 days 
Fine Grid) At= 9 hours 
Atmosphere Only 697 days 
(Fine Grid) 
Combined Ax=Ay=200km 140 days 
(Linear) At tm 30 min 
(Not used) AC = 9 hours 
ocn 
linear momentum 
equations 
Combined Ax=200km 320 days 
(Mixed Grid) Ay=400km 
(Not used) AE =] hour 
atm 
At =9 hours 
ocn 


TABLE 2. Review of computing time requirements for 
variations of the numerical models used in this experiment. 
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It can be readily deduced from these time figures that 
Me coupling processes required lengthy calculations. The 
method by which the sea surface temperatures are transferred 
Into wave number space required many time consuming opera- 
tions for every oceanic time step. The other summations 
and averaging procedures all contributed to a joining pro- 
cess that caused the combined model to run significantly 


slower than either of its individual components. 
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VI. ZBIESEUSSTON 


The data that was collected during the long term inte- 
gration can be placed in two categories. The full fields 
of data and the data collected at the five designated 
points differ widely in frequency of observations and area 
of coverage. Because of this, the methods by which the 
data were to be examined differ considerably. The full 
fields of data that were retained at 100 day intervals were 
examined in a synoptic-type study, while the five-point 
data, recorded at daily intervals, were examined using 
spectral techniques. 

In the synoptic study the analyzed data were compared 
not only to observed climatological data, but also to the 
results obtained from the various phases of this numerical 
experiment. The atmospheric energy balance was examined 
and comparisons were made between phase I and phase III and 
with the work discussed in Chapter II. In the spectral 
Study,power spectrum estimates that were obtained from data 
at one particular point were compared with estimates ob- 
tained from data at other points during the same phase of 
the experiment. Similarly, spectral estimates of one 
variable were compared with estimates of the same variable 
at the same point in a different phase of the work. Some 
comparisons were made with spectral estimates calculated 
from observed data collected by the Woods Hole Oceanographic 


Institution. 
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BOSSE SYNOPTIC STUDIES 
Throughout the experiment full fields of the stream func- 
won, vertical shear current components, and temperatures 
were retained as oceanic data and the Fourier coefficients 
were saved as the atmospheric data. These variables, along 
with fields of oceanic vertical motion and atmospheric zonal 
velocity, vertical motion, mean meridional velocity and 
disturbance amplitude and phase information were printed 
every 50 or 100 days. Many of these variables were examined 
in detail and discussion of the results is contained in the 
subsequent sections. 
1. The 400 km Grid Combined Model Experiment 

The 400 km grid combined model was merely a tool that 
was used to expeditiously integrate the model equations to 
reach a state of quasi-equilibrium. The ocean-atmosphere 
interaction at the end of this part of phase I was not of 
primary interest. The state of the system at the end of 
day 20,200 (55 years) was to be used only to initialize the 
fine grid model from which more realistic and more meaning- 
Zulzresults may be obtained. However, a brief discussion 
Pssehnescomdıitions to which this form of the combined model 
evolved is necessary to explain why the 200 km grid was 
needed to complete this phase. 

Micwtnrietal -condLelons. for this parteof the “esperas 
cie mp evilously been described. The initial adjustment 


of sea surface temperatures was due to the large horizontal 
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May drrtuslom epesficeient then in use. The redwetion of 
this coefficient and the necessity for maintenance of a 
large meridional sea-surface temperature gradient was al- 
tady discussed. The coefficient of lateral diffusion 
could not be made small enough with the 400 km grid so that 
nonlinear terms would be important. Bryan and Cox (1967) 
Mave detined an “effective” Reynolds number, Re, as the 
Reynolds number associated with the lateral turbulent mix- 
ing. The nonlinearity of the solution is dependent on the 


size of the Reynolds number defined as follows: 


U*L 


> 
Ay 


Re = 





(619 


where U* is the velocity component depending on the density 
gradient, 
L is a scale length, and 
Ay tthe Coefficient of lateral momentum diffusion: 


Bryan and Cox conducted a series of experiments 


with values of Re = 0,8, which gave a highly diffusive, 
linear solution and Re = 8.0, which yielded a moderately 
diffuse, slightly nonlinear solution. The maximum value 


of Re used in this experiment was 6.5 and it is believed 
that a value greater than 10 was necessary to achieve a 
truly nonlinear solution. It is easily shown that with the 
AMO @raid, continued reduction of Ay produces a situation 
where the linear western boundary current becomes unresolv- 


mote poee, Equation (4.1)] prior to the occurrence of non= 


linear effects. Therefore, the solution with the 400 km 
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erid is similar to the viscous linear solution described 
PrsMunk (1950). Figure 21 shows the ocean stream function 
ebtained at day 20,200 (about 55 years). According to 
Munk, the streamlines should represent an exponentially 
Mecavine Series of vorticies centered on the latitudinal 
axis. The vorticies are present across the width of the 
domain but do not appear to decay as expected. Truncation 
error is a probable cause for insufficient decay in this 
area. In the gyre north of 45N, the maximum value of the 
vertically averaged meridional current was on the order of 
TO cm een. In the southern gyre the maximum value of the 
Same quantity was about 15 cm see. The magnitude of the 
Beraleskwertical shear plus vertical mean parts) meridional 
currents in the northern and southern gyres were 10 and 35 
cm sec respectively. The higher latitudes in the model 
were more barotropic and the vertical mean current domina- 
ted, while the vertical shear current dominated in the 
Rc baroclinic areas of the central portion. In the 
Borcehschesvertical shear current was about zero, while in 
the south the vertical shear current was about 20 cm sec 
directed northward. 

Figure 22 shows the sea surface temperature field 
ALI The mean meridional temperature gradient 
over the latitudinal domain was 21.4 degrees per 6000 km. 
The temperatures at both the north and south boundaries had 
decreased from their initial values. Upwelling along the 


Pome naeeportion of the northern boundary was responsible 
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Latitude 


18% 


Longitude 


Pearce Ocean stream function obtained 
from 400 km grid combined model. 
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at day 


20,200 





72 


Latitude 


181 


Figure 22, 


Longitude 


Sea surface temperature at 
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day 20,200, 
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Acs cooling at the surface. The small scale spatial 
fluctuation in the temperature field in low latitudes ap- 
meatreewto be a result of the linear viscous solution in 
which there is alternating northerly and southerly flow and 
therefore northerly and southerly advection of heat. 

Figure 23 shows the zonal winds at three levels in 
the atmosphere. The curve labeled (a) depicts the zonal 
mean surface wind. The easterly, westerly, easterly pattern 
from north to south is qualitatively like the patterns in 
the real atmosphere. However, the latitudes at which the 
mercetion of flow change occur are slightly north of the 
actual latitude of change. Curve (b) is the zonal wind at 
750mb and curve (c) is that quantity at the 250 mb level. 
The increase of westerlies with height in the central 
latitudes is in thermal wind balance with a temperature 
Serwetuüre of warm air in the south and cold air in the 
North. 

While the structure of the atmosphere-ocean system 
Mat evolved from integrating the model equations over the 
so km grid was not extremely accurate, it did suffice for 
initializing the integration of the fine grid model. 

e Om Grid. Combined Model Experiment 

rhe fine grid or 200 km grid part of phase I was 
eaıtsalsrzed with data from day 20,200 of the 400 km grid 
mens Bersspatjalrsinterpelation. The model equations 
were then integrated in time for 8700 days (approximately 


24 years). The adjustment to the fine grid was evident in 
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a slight initial decrease of the oceans volume averaged ten- 
perature. After a few hundred days, this temperature start- 
Pemiunecreacing and at day 21,000 the sea was warming at a 
fate of 0.8 degrees per century. This warming rate was main- 
tained for over 1000 days and then the rate decreased grad- 
Have to 0.5 degrees per century at day 28,900, the last day 
of the experiment. This warming rate was caused by a down- 
Mama heat flux across the air-ocean interface of about 2.3 
cal aa, While the experiment could have been con- 
tinued until no heat flux was observed, it was believed that 
the large amount of computing time required to accomplish 
this would not have been worth the slight change in the 
final results. 

At the end of the 400 km grid model experiment, the 
total disturbance energy in the atmosphere was statistical- 
ly steady, oscillating about 50 ae SN During the first 
few e hundred days of the fine grid integrations the average 
disturbance energy jumped to about 70 me ec but then 
quickly returned to the statistically steady state which 
kad previously been observed. The initial cooling of the 
ocean imparted an added amount of heat to the atmosphere 
in the southern half of the region generating more zonal 
Potential energy which the baroclinic processes converted 
into disturbance energy. After the cooling ended, the 
Warming of the ocean apparently occurred rather uniformly 


in space resulting in the generation of about the same 


amount of available potential energy as had been generated 
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waon ieo reducing the grid size. More discussion of the 
Eme ey balance in the fine grid model is in a later section. 

The mean state of the atmosphere over the last 100 
days of this phase of a Eonia 15 shown in Figures 24 
Ao These maps along with Figures 28 to 30 that will be 
discussed shortly are deviation values of the height fields 
Various levels obtained from the geostrophic condition. 
Mie zero contour and plus and minus deviation contours are 
drawn. These maps are displayed in this way because the 
@eviation of the stream function is predicted by the quasi- 
geostrophic atmospheric equations. The quantity Y appears 
2U in all terms, and can be represented as V(Y'"(x,y,t) 

T V(t)), a deviation and mean part of the stream function. 
Therefore, any mean could be added to (or subtracted 
mom) the deviation quantity. 

Figure 24 shows the extrapolated distribution of 
the 1000 mb height deviation analyzed at 12-meter intervals. 
There is a semi-permanent mid-latitude cyclone over the 
eastern ocean and a high at the same longitude south of 
the low. A high is present in the northernmost section of 
Eres region, This is worthy of note since there is no cold 
land mass modeled into the system to simulate a "Greenland" 
wne effect. Figure 25 shows the distribution of the 250 
mb height deviation analyzed at 60-meter intervals. The 
Mes zonal throughout the central portion of the domain 
Paieaeey;clLonic polar vortex present in the north and a 


weak ridge present at the same longitude in the south. 
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Figure 24. Distribution of 100 day mean 1000-mb height 


deviation in phase I. 


AZ 





72 


Latitude 





Longitude 


Pureza. Distribution of 100 day mean 250 mb height 
deviation in phase l. 
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Be 26 shows the distribution of the 750 mb height de- 
pato analyzed at 30-meter intervals. The flow is again 
zonal throughout the central portion of the region. How- 
ever, the lowest height is now in the northeastern corner 
of the area and the low latitude high is present in only 
ebezssouthernmost part of the domain. 

The vertical structure of the standing disturbance 
is such that the 1000 mb high in the polar regions lies 
below the 250 mb low the same as in the real atmosphere. 
The warm low pressure area at the southern boundary at the 
1000 mb level becomes an area of high pressure at the 250 
mb level. 

The zonal wind averaged over the same 100 day per- 
iod is shown in Figure 27. The line labeled (c) shows a 
westerly flow at all latitudes at the 250 mb level with the 
ietmbocated at 45N, coinciding with the strong height 
gradient in that region as seen in Figure 24. Curve (b) 
shows the weak polar easterlies at /50 mb level to the 
north of the low at 60N. The remainder of the flow is 
westerly with a weak jet at about 45N. Line (a) shows the 
zonal flow at the 1000 mb level. The polar easterlies 
extend from the northern boundary to a latitude just north 
Wim mec rontece activity  Weèsterlies extend from that 
Pome to the latitude of the 1000 mb highs and easterlies 
cone southern boundary. <A comparison of Figure 27 
with Figure 23 shows that the mean zonal winds did not 


eenei geniiicantly after reducing the grid size. This 
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treo Distribution of 100 day mean 750 mb height 
n ation in phase I. 
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N entes that for the atmospheric model, a 400 km grid 
would have sufficed with little loss of resolution. 

Figure 28 shows the distribution of the 1000 mb 
height contours at 25 meter intervals at day 28,900. A low 
dominates the western half of the region at this time and a 


high covers the eastern half. This figure was compared to 
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Pure 27. Zonal atmospheric winds in m E from 200 km 
grid model at a) 1000 mb; b) 750 mb; c) 250 mb. 

the series of 1000 mb figures presented by we hilltop sae 95169 
and the gross features are almost identical to the early 
stages of his model atmospheric development. The main sur- 


face troughs and ridges are orientated so as to lean back 
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remazcethe northwest in the northern half and slightly to- 
ward the southwest in the southern half of the region. The 
eomversgence of momentum in the central portion increases 
the westerlies and accounts for the mid-latitude jet ob- 
served in this area (Figure 29). 

paez shows the distribution of the 250 mb 
height contours at 60 meter intervals for day 28,900 (79 
years). A wave of moderate amplitude is superimposed on 
the zonal flow with a trough present over the western half 
of the region and a ridge dominating the eastern half. The 
polar vortex observed in Figure 25 exists at this time 
also. The single wave structure of the atmospheric model 
merecesma pattern similar to this to be found at any given 
instant. The disturbance moves and its amplitude varies 
with time from a zonal pattern to a highly amplified wave 
Battern. At this time the disturbance kinetic energy had a 
value of 54 me ee This is quite large compared to the 
disturbance kinetic energy of the more zonal flow in Figure 
25. However, at other times in the experiment higher 
values of this quantity were observed. It is believed that 
the atmospheric state shown in Figure 28 is representative 
of a flow in which a wave disturbance of average amplitude 
momdimbedded, At times the amplitude of the wave will be 
Seater than it is at day 28,900 and at other times less. 
Also, the wave will propagate across the domain and yield 


a time mean upper level flow as shown in Figure 25. 
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Figure 28. Distribution of 1000 mb height deviation at day 
Mo 00 in phase I. 
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ez Distribution of 220 mb height deviation at day 
900 in phase I1. 
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ceo Ene distribution of the 750 mb height 
eont our at JUemeter intervals at day 28,900, shows a closed 
Sec lommcneclrenlation directly beneath the trough at 250 mb, 
The amplitude of the wave motion is about the same as at 
250 mb and a closed high is present beneath the upper level 
ridge. 

The suggestion that the single wave in the x-direc- 
tion would produce a series of cyclonic scale disturbances 
moving over the domain of the atmosphere-ocean model was 
born out by these figures plus examination of preceding 
maps. The wave moves eastward with a speed of about 600 km 
per day. This is considerably slower than the wave motion 
echallıips' experiments. 

The vertical heat flux was discussed earlier as be- 
ing proportional to the atmosphere-ocean temperature differ- 
ence. However, due to the method by which the sensible 
heat exchange is calculated, the vertical heat flux is 
actually proportional to the difference between the atmos- 
pheric surface temperature and the Fourier analyzed sea 
surface temperature. This temperature difference, averaged 
Bere lasıt 100 days of the experiment, 775 shown in 
Figure 31. In general, there is cooling of the atmosphere 
in the western sector and warming of the atmosphere in the 
ac a le this heating may be representative of what 
eremrs ane the northernmost latitudes of the real atmosphere, 


does aL accurately portray what occurs throughout the 
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Figure 30. Distribution of 750 mb heicehtedevsatrionratsd., 
28,900 in phase I. 
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Figure 31. 100 day average difference between the Fourier 
transformed sea surface temperature and the atmosphere sur- 
face temperature computed in phase I. Heating of the 
Semosphere Occurs where values are positive. 
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peor wey OL the Oceanic régions where cold atmospheric air 
masses move from over a cold continent to a position over 
Naci water and are heated from below. The lack of colder 
continents in the model and the use of cyclic boundary con- 
ditions makes it impossible to simulate this phenomena in 
this simple model. The use of the Fourier analyzed sea 
Surface temperature also contributes to the unrealistic 
vertical heat flux pattern of the model. The strong tempera- 
ture gradient observed in the real ocean along the east 
coast of continents and observed in the model along the 
western boundary of the ocean in Figure 34, had little 
effect on the heat flux pattern that evolved. The Fourier 
analysis smoothed this gradient and although the gradient 
did influence the new sea surface temperature somewhat, it 
did not have the same influence observed in the real atmos- 
phere-ocean system. The Fourier analyzed sea surface tem- 
perature was used in the heat exchange to be consistent 
with the simple one-wave formulation of the atmosphere. 

memc yclic boundary condition forces the air to pass from 

a region of warm water on the eastern boundary immediately 
to the region of cooler water in the western area. The 
overall effect that the resulting heat flux pattern had on 
Aita obs of this experiment are hard to determine, Ine 
strong heating observed in nature on the western side of the 
ocean in the mid-latitudes does not occur. The overall 
Sriects may be to smooth out the effect of the heating on 


Mew rtimospicric temperature disturbance (only disturbance 
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quantities are of interest since domain mean temperatures 
Orsehezatmosphere do not change) and cause a more zonal 
partern. Figure 32 shows the time mean atmospheric surface 
temperature averaged over the last 100 days of the combined 
model experiment. The temperature in this figure has a 
more zonal appearance than was expected after examining the 
meaneoO0 mb height deviations (Figure 24). 

Figure 33 is the 100 day time mean of the stream 
Eunetilon governing the vertical mean currents. The sub- 
arctic gyre, visible at about 55N in this figure, produces 
southward flow along the western boundary at that latitude 
and northward return flow 4-600 km east of the boundary. 
Mets ab tropical circulation at about 35N has northward flow 
at the western boundary and southward flow 4-600 km to the 
east. The convergent flow in the boundary causes a sharp 
meridional temperature gradient at about 43N and contri- 
butes to the low temperature near 48N in the western bound- 
ary. 

Figures 34, 35, and 36 show the 100 day time means 
of the ocean temperature at the surface, 500 meters and 
1300 meters, respectively. These maps are vertically con- 
sistent and the low temperatures just described near 48N can 
be seen at all three levels. The thermal ridge observed 
ime tiem monrthiwestern sector of Figure 34, and also in Figure 
35, is caused by the northwestward flow on the east side 


of the subarctic gyre and somewhat by the upwelling along 
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pare 34. 100 day mean sea surface temperature (C) 
paterno in phase LI. 
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Figure 35. Temperature pattern at a depth of 500 meters in 
aean in phase I (C). 
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Figure 36. Temperature pattern at depth of 1300 meters in 


Pan phase 1(C). 
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miemwestern boundary seen in Figure 37, which contains the 
ame averaged vertical motion field at the depth of 300 
meters, In the fine grid, as well as in the coarse grid 
discussed earlier, the flow in the western boundary regions 
in high latitudes is composed almost entirely of the verti- 
ea mean current. Here the flow is of this character and 
Moves across the isotherms, advecting warmer water toward 
the northwest. The large area with horizontally isothermal 
meter Of about 5C located west of the ridging in the north- 
west quadrant of the ocean domain can be compared qualita- 
tively to a similar feature found east of Newfoundland in 
the North Atlantic. Sverdrup et al (1942) presented sea 
surface temperature maps that contained this feature in 
both winter and summer. 

In the northeastern quadrant the flow is weak and 
aio se ly parallels the isotherms. Here, in the absence of 
empreciable vertical mean current, the vertical shear cur- 
rent is dominant and when the flow reaches the northeastern 
boundaries large sinking motions occur (see Figure 37) pro- 
ducing warm temperature at all depths there. Figures 34, 
Soeand 36 ali show warmer water in the northeast corner 
due to this phenomenon. The isotherm pattern produced by 
the model in this area can be likened to the real sea sur- 
Acre nperature pattexn off the west coast of the United 
States. The area of warmest surface water in Figure 34 is 
observed in the southwestern corner of the domain. This 
aredd is analogous to the Sargasso Sea with its warm tem- 


peratures and weak currents. 
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Figure 37. 1007 day’ mean vertical motions in cm day 
300 meters in phase I. 
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The temperatures along the western boundary are fur- 
ther affected by the upwelling that takes place due to the 
Bermalseirenlation. The central latitude is a region of 
weak eastward flow in the surface layer and westward flow 
at the lower layers. The moderate upwelling brings colder 
water to the surface, lowering the surface temperature 
along the boundary. An area of weak sinking motion in the 
central portion of the domain in Figure 37 brings 
about the warm 8C water at the 500 meter level (Figure 35). 

Paemetie:stream function field in Figure 33, the 
merttcal mean current can be computed. The subarctic gyre 
has a vertical mean current on the order of 20 cm sec 
The subtropical gyre has a vertical mean current on the 
order of 15 cn ee Table 3 shows the total meridional 
Peolocity, an estimate of the vertical mean meridional cur- 
Sent as calculated from the stream function and the verti- 
Gal shear meridional current calculated as a residual, in 
“be western boundary current. The vertical mean current 
is almost entirely responsible for the meridional current 
the boundary layer associated with the subarctic gyre, 
poten the vertical shear current is dominant in the sub- 
ErOpical gyre. The maximum value of 51 cm ec in the 
Porimeaeyseurrent, repnresenfing a Gulf stream type flow, 1s 
less than that found in the real world oceans where 100 
cm sec” üevesbeen observed for this type Ons a boundary 


Current. 
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Latitude Total Merid. Vertical 
Current Mean 


Vertical 





Shear 
2.60 ZO -0.2 
12.8 O a 
>, OO -0.1 
-0.4 -1.7 S 
-12.0 -14.0 20 
-20,5 -19.0 -1.5 
-18.7 -14.0 -4.7 
-5.6 26.3 0.9 
EL 0.8 7043 
E7 Sc TOTS 
SD eS, 250 
Deel 5:75 Sn 
21,0 8.0 33.0 
2055 239 1920 
6.7 =5. 1l Lrg 
w9 -0.7 2:6 
DABLE 3. Total meridional current (cm a as produced 


by the combined model; vertical mean current as computed 
won the model stream function; vertical shear current 
calculated as a residual. 
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Figure 38. 100 day mean stress pattern produced by con- 


bined model. 


mestress acting as awdr yimestorce son the ocean 
vas Continuously varying as the time integrations progressed. 
A 100 day time mean zonal stress was calculated at a time 
six years prior to the end of the integrations and is shown 
in Figure 38. This stress was used as a constant force 
eonte ocean during the ocean-only phase of the 
pe me NOW it can be anterpreted as the mewn stress 


esponsble for the stream function field observed in 
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Beer. une curl of the stress has maxima at about 34N 
and 50N which are the approximate positions of the centers 
of the subtropical and subarctic gyres,respectively. The 
positions of the maxima are about 5 degrees further north 
than actually observed in nature. This is responsible for 
the northward shift of the gyres and their associated 
boundary currents. The stress varies from 1.20 to -1.15 
dynes aa These are reasonable values lending justifica- 
tion to the stress formulation discussed earlier. 

The state to which the combined model progressed 
was acceptable in its resemblance to the real atmosphere- 
ocean system. While it was not the intention of this work 
to develop an operational model for forecasting either 
atmospheric or oceanic parameters, it is comforting to See 
a reasonable distribution of the varlous parameters. 

3. The Ocean-Only Experiment (Phase IT) 

The zonal mean stress (Figure 38) that had been 
averaced over the 100 day period from day 26,650 to 26,750 
Wassused as a constant forcing function for the subsequent 
ei ear period over which the ocean-only model equations 
merecen teserated. The atmospheric Surface temperature, Ty» 
was also averaged over this same 100 day period. Figure 
Ros ta plot of the zonal mean temperature as a function 
of latitude. The temperature varied from 2./ degrees 
dental rade in the north to 23.3 degrees in the south. This 
deneraat üre was used as an equivalent atmospheric temperar 


ture to which the sea surface temperature was compared to 
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determine the vertical heat flux. Except for the fixed 
Memosipieri.c surface temperature which was independent of x, 
the same formulation was used for the heating as was used 
in the combined model. Now, however, one degree of freedom 
was removed and the ocean in avery short time became 


paeeady, except for a slow secular change. 
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w e anelon of Latitude. 


The instantaneous fields of the ocean dependent 


Variables from day 26,700 were used as initial conditions 
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for the ocean-only phase of the experiment. During the six- 
year ocean-~only experiment, data fields for all oceanic 
variables were saved every 100 days and selected daily para- 
meters were saved at the same five points as previously 
discussed. 

The 2200 day (six year) integration of the ocean- 
only equations proceeded at a rapid pace. No problems were 
encountered with adjustment to the new constant forcing 
functions and a quasi-steady state was achieved. At the 
Start of the 2200 day phase of the experiment the volume 
average temperature of the ocean was increasing at about 
0.5 degrees per century as previously discussed. Through- 
out the six year period the warming rate remained about the 
Same, with a 0.5 degree per century warming still present 
at the end. This makes this author believe that perhaps 
aS S as close to an equilibrium state as could be 
reached with the model as formulated, without the expendi- 
bWweemar an excessive amount of computing time. An examina- 
tion of the sea surface temperature, Figure 40, and the 
area average of the temperature at that level shows that 
the sea surface temperature is not changing at all. The 
rate of change of average temperature at the upper level 
EE e tly zero. However, the average temperature arte 
ovest level is changing at a rate of 0.6 degrees per 
Century, showing that the diffusive mechanism by Which Ehe 


deeper water must be changed is very slow. 
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Pure 40 shows the sea surface isotherms at day 
26,900. An examination of this temperature field at day 
28,800 was made and the two maps were identical. Therefore, 
it was considered unnecessary to derive time mean maps of 
Piieoedwantity. The pattern is similar to the sea surface 
temperature field at the end of phase I (Figure 34) in that 
a strong temperature gradient is present in the western 
boundary layer at about 35N. The isotherms throughout the 
emer ale portion of the domain are slightly more zonal titan 
in phase I where there is an x dependence in the forcing, 
Dmemtie meridional gradient in this area remained about the 
same. “The warm core tropical water (above 23C) does not 
extend into the boundary layer but is visible as a narrow 
ima ot warm water just north of the southern boundary. 

Ohe major differences instheytemperarcune, Liclds ves 
tween the two phases is in the northern half of the domain. 
The thermal ridge present in the northwest quadrant in 
phase I is observed only as a small perturbation in the 6C 
isotherm in this phase. The thermal trough to the west of 
As eature is similar to that present in the combined 
model with a region of water with temperature between 4 and 
6 degrees extending over a thousand kilometers. The remain- 
der of the isotherms in the north have a more zonal appear- 
wee here is a small pocket of water below 2C in the 
northernmost region in phase LI that is not present in 
phase I. The isotherms in the northeast sector vary signi- 


ficantly between the two phases. In phase I the thermal 
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rent 2 lows into the northeast corner when the strongest 
remem motion occurs. , In phase II the vertical shear cur- 
rent flows into the eastern wall on a more eastward heading. 
The maximum sinking motion is located about 800 to 1000 km 
Sot its position in phase I. This can be verified by 
eompating Figure 41, the vertical velocity at a depth of 
300 meters in the ocean-only case, with Figure 37, the same 
quantity in the combined model. The northward curvature of 
the isotherms in the central portions of the eastern bound- 
ary 1s present in both phases, but is less pronounced in 
the ocean-only case. An examination of the meridional cur- 
rent in this area showed that there was weak northward 
motion in the northeastern boundary region in both phases 
Which when coupled with the downwelling in that area, was 
responsible for the warm water there. The temperature 
patterns for the levels at 500 and 1300 meters resemble 
their counterparts from phase I (Fiure and 30) vithe 
same basic variations as were seen in the upper level. The 
maximum in the temperature in the northeast corner has moved 
south about 800km as the region of sinking motion moved. 
The isotherms show no unexpected results and will not be 
depicted here. 

Miér vertical motion patterns play an important pare 
Meethe wad ;ustment of temperature patterns in the boundary 
Holmes Erong vertical motions occur there due to the 


thermal circulation. A region of sinking water near 40N 
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in the western boundary region causes the weak gradient be- 
tween the 12C and 14C isotherms. The isotherms along the 
western boundary are shifted southward. That is, the 14C 
isotherm in the ocean-only model is at the latitude at 
which the intermediate 15C line (not drawn) was in the com- 
bined model. 

Figure 42 shows the stream function field for the 
ocean-only model at day 28,900. It was also found unneces- 
Sariy to compute a time mean of this quantity. The subarctic 
eyre and the subtropical gyre both appear to be slightly 
further north in the ocean-only model than they are in the 
combined model (Figure 33). Those gyres have about the 
same intensity in the western boundary region as the corres- 
ponding gyres in phase I. Tnis was expected since a mean 
stress from the combined model was their driving mechanism. 
The northwestward flow on the east side of the subarctic 
gyre is still present and causes some warm water to be 
advected, but at this time the large amplitude ridge that 
was observed in phase I was not present. The ocean-only 
sea surface temperature appears to be strongly tied to the 
imposed zonal atmospheric temperature, except in the west- 
Sabio in dary region and to a lesser degree in the eastern 
Boundary region. 

Ehe maximum value of the meridional current 1s re- 
Miccamoonewiat in the ocean-only model, In the western 
boundary region, the maximum meridional velocity is 47 cm 


= a. 
sec G in the ocean-only model compared to 31 cm sec in 
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ena. “Stream function pattern at day 28,900 in 
phase II. 
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the combined model. Since the vertical mean current appears 
dment tehe same in both models, the difference is in the 
aeea near current. The difference in the predicted 
ocean variables between the ocean-only and combined model 
phases is small. Radical changes did not occur. The dif- 
ferences that are present are due to the removal of the 

time varying atmospheric forcing mechanism in both the ex- 
change of momentum and heat in the ocean-only model. 


4. The Atmosphere-Only Model (Phase III) 


Mme zonal part, a of the combined model sea- 
surface temperature was averaged over the same 100 day 
period as were the zonal stress and atmospheric surface 
temperature, This quantity was used as the underlying tem- 


perature with which the atmospheric temperature was com- 


Baredgsfor the calculation of the vertical flux of sensible 





heat, Figure 43 shows the profile of the time and longi- 
ie 
midiial mean of sea surface temperature, T * which 


sea 


closely resembles the mean atmospheric temperature shown 
AFi gure 39. 

imemwanitial conditions for phase Lil were the 
Potter coetficrents, A through F, from day 26,/00 of 
phase I. There was essentially no adjustment required and 
PBuer2200 day integration was conducted, with data retained 
as before, in a minimum of computing time. At the end of 
this phase the mean state of the atmosphere over the last 
100 days was examined along with the state of the atmosphere 
From day 28,900. These maps were compared with the atmos- 


phere at the end of phase i. 
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Pagure 43. 100 day mean zonal component of sea surface 
memperature (C). 

Figure 44 shows the 100 day mean distribution of the 
MODO ab height deviation as a function of latitude. Since 
iemenermal forcing is now a function of y only, the mean 
fields are purely zonal (with the exception of small errors 
due to averaging and truncation). The deviation in height 
extends from 18 meters at 35N to -28 meters at 54N and the 
deyiarıon pattern is Such that a geostrophic wind regime 
of easterlies, westerlies, easterlies would exist (as it 


does). Figure 45 shows the 100 day mean distribution of 
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Figure 44, Distribution of 100 day mean 1000 mb height 
deviation as a function of latitude. 

the 750 mb height deviation as a function of latitude. The 
deviation goes from 118 meters at the southern boundary to 
-120 meters at the northern boundary. The tight negative 
height gradient through the central latitudes agrees with 
the existence of a regime of westerlies in that area. 
Figure 46 depicts the distribution of the 100 day mean 

Ae tr ibution of the 250 mb height deviation as a function 


Semplatitude. The gradient that exists from 410 meters in 
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er zsouch to-407 meters in the north also agrees with the 


westerly regime at that level. 
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mene 45, Distribution of 100 day mean 750 mb height 
den ation as a function of latitude. 

Instantaneous maps at day 28,900 of the distribution 
e neight deviation show a highly disturbed atmosphere: 
eeure; the 1000 mb height deviation field at day 28,900 
analyzed at 25-meter intervals, shows a pattern that has 
a disturbance amplitude which is about the same as its 
counterpart in the combined model (Figure 28). Figure 48 


Smowsmthne distribution of the 250 mb height deviation 
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analyzed at 60 meter interyals. This flow is quite similar 
mehat at the same leyel at the same time from the combined 
modell {Pigure 29). Figure 49 depicts the distribution of 
Bue 730 mb height deviation analyzed at 30 meter intervals. 
The wave-like structure here has a greater amplitude and 

Bes deviation is larger than the same fields in the combined 
medie (Figure 30). The occurrence of the difference in the 
mipettude of the patterns is not surprising since the model 
memosphnere is continuously changing. It would be a coin- 
peiidenece if the structure of both the atmosphere-only and 


combined model atmosphere were identical at any given time. 
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Eee Distribution of 100 day mean 250 mb height 
ataron tas a function of latitude. 
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The 100 day average zonal wind profile in the atmos- 
phere-only model is very similar to that in the combined 
models hmormeinskRieure 27:2 Avery slight reduction in the 
jet maximum at 250 mb is the only difference so the atmos- 
Piore—on | yevimcdeprott te willenet be presented graphically. 

The time averaged results of the atmosphere-only 
phase of the experiment were very similar to the time 
averaged atmospheric state arrived at by the combined model. 
Since the sea surface temperature is a slowly varying para- 
meter compared to the atmospheric time scale, holding it 
constant approximates long term actual conditions fairly 
A EE ro ula be surprising if a radical change occurred 
in the time averaged fields in the atmosphere-only formu- 
Eat ion: 

one inn yeaa 

An examination of the energy in the atmosphere was 
made at the completion of phase I and phase lll. Rather 
than study instantaneous values of the energies and energy 
transformation as was done in the earlier discussion of the 
energy balance, values averaged over the last 100 days of 
Ene experiment were examined. The energy flow convention 
as established in Figure 4 was used. 

The mean energy exchanges for the combined model 
and the atmosphere-only model are shown in Figures 50 and 
Dine respectively. EE eet erestinge to note the similarity 
in the flow diagrams in the two phases. The diabatic heat- 


ing generates about the same amount of zonal potential 
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Hure 50. Plow diagram of mean atmospheric energy from 
the combined model. Energy and ae transformation 
Teer nva m ( = 10 "ere cm-- sec’). 

Ven =r he bOLmmeaces . lwo Watt Ho mie pilacew ss and g 
watt m Picci te lnemaonversionnor zonal potential 
me vweOoce Mt alsendschensinto eddy kinetic energy is 


: : : -2 
M e a n boti phases at a value of 1,26 watt m À The 


Conversion of eddy kinetic to zonal kinetic energy is also 
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i i 
Figure 51. Flow diagram of mean atmospheric energy from 
the atmospheric-only model, Same units as in Figure 50. 
very similar in the two phases with values of 0.37 and 0.33 
watt nes tiemdiccipatiensomethe Kinetic energy (both 
aona ae day by Boch internal and Surface friction are 

IA co lla wa ie wOnbyeconvyersion that differs by a 
significant percentage is the conversion of zonal kinetic 


to zonal potential energy. However, this has been pointed 
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out earlier as being a very small quantity and no comment 
on the difference between the two phases is offered here. 

It is not surprising that these flow diagrams re- 
semble each other. The only difference in the two models 
lies in the method by which heat is exchanged between the 
atmosphere and ocean or underlying wet surface as the case 
(iwc. En the calculation of the sensible heat flux 
necessary for the generation of zonal potential energy in 
the combined model only the zonal part of the sea surface 
temperature was used. Therefore, when the disturbance 
part of the sea surface temperature was neglected in the 
atztmosphere-only model, it did not affect the generation of 
energy very much. In the other energy exchanges the differ- 
ences in the disturbance quantities were not sufficient to 
make much difference, 

These diagrams compare favorably with Oort's dia- 
gram (Figure 6) as did the earlier energy diagram (Figure 
Dee lhewdissipatiton of “eddy kinetic energy is still too low 
wo new explanation is offered for this. It is somewhat 
Surprising that the instantaneous values of the energy ex- 
change shown in Figure 5 compare so favorably with Figures 
Jia so many vays. The fact that the time rate of 
ehane or kinetic and potential energy (both zonal and 
eddy) were other than zero in that case shows that either 
there were small errors in calculations, or an exact steady 


SlLAabcet iwethemmencr ea ies did not exist, as was stated in the 


199 





Barilıer discussıon.  Th2s is in accord with present know- 
E o tne neal atmospmere which is never in exact steady 


State: 


BSZESEECTERUM ANALYSES 

This section will briefly discuss the spectrum analysis 
Beehniques used in this study. A computer program in the 
BROMEDTICAL serres (BEMED O21) developed at the UCLA Medical 
Genter (Dixon 1970) was modified and utilized for all 
spectral analyses conducted on the data generated in this 
study. 

Basic Derinitions 

Prior to discussing the procedures used or the 
results obtained, a brief review of the definition of a 
einst he basic statistical parameters is given. 

a. Variance - the mean square deviation from the 
mean of a variate. 

bem tcandard deviation -~ the positive square root 
of the variance, used to measure the degree of variability 
ol a variate. 

O u eo Correlati on coeficients ~ linear correla- 
tion coefficients between a time series and the same time 
series an interyal of time later. 

to covariance ean auto-correlation function 
that has not been normalized by the yariance. 

E oS Correlation - correlation between two 


mece es Monero nen hs at an interval of time later. 
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Tesserose ceovarıance = a cross=-ceorrelation function 
Eharshas not been normalized by the variance. 

AO ES Decano time Series - a function of 
Ereomerevouwhichspivwessthezseontribution of oscillations at 
facie requencyerorthezsvsrıance of a time series. 

Domsstcormeekerume- S2srunetion of frequency which 
nes the contribution Of oscillations at that frequency 
tothe total cross-covaáariance of two time series. 

II es peana eunetion of frequency 
Cech pei ves the contribution of oscillations at that fre- 
queney to the total covariance of two time series X(t) and 
Wt) obtained when all the harmonics of time series X were 
delayed by a quarter period but the Y series remained un- 
changed. 

ae tne, relative phase of the harmonics 
and Y(E) is measured by the ratio of the quadrature 
meetrum to the cospectrum. 

2% Spectral lechiniques 

Many investigators (for example, Yanai et al, 1968; 
and Wallace 1971) have utilized spectral techniques to 
study tropospheric phenomena primarily in areas of sparse 
data where a single station's time varying record is 
Pa ne melszeiınportant.sllersamestype of spectral techniques 
Came berusedero Study thesatmospheric and oceanic data ob- 
Coe ot he numerical integrations conducted in this 


Sua) MoO e ive selected points (Figure 20) six years of 
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Beer eceorlded ar 1 izo day intervals. At each point 
1960 observations were obtained. These data were analyzed 
and a power spectrum was obtained for most variables at all 
five points. Cross-spectrum analyses were conducted on sevyv- 
eral variables yielding cospectra, quadrature spectra and 
phase relationships. 

More attention was focused on point C because of 
its position in the western boundary region, It was located 
in an area where meridional velocities were strong, tempera- 
ture gradients were tight and upwelling was important. 
However, the other points in Figure 20 were not neglected 
and comparisons of the resulting spectral statistics between 
points were made. 

The data were recorded at an interval, (A t) a> of 
1.125 days (three oceanic time steps) during the three six- 
year experiments. That is, the data were available for 
analysis from the combined, ocean-only and atmosphere-only 
pmases Of —thewvexperiment. the decision to record data at 
w days was a critical one. It would require an exces- 
Sive amount of storage space to keep a longer record of 
daea at Shorter recordinerinteryals and problems could be 
eweountered ir the interval was too long. It is important 
to ensure that the series which is being examined does not 
contain substantial amounts of variance from frequencies 
outside the range that can be resolyed with the spectra, 
spectral information for frequencies less than one cycle/ 


2(Lt),; the Nyquist folding frequency, can be obtained from 
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nel re eouencireschicher than this are “folded” or 
Ia stitomritcm lower trequency region. Wallace (1971) 
Ders sous gez recneshrehutrregqueney 131mit does not appear to 
pose any serious problem from aliasing because there is 
Member varuancesassocrated with frequencies higher than 0.5 
Salas per dayveinwtmemamuermediate to large scale disturb- 
ances. From experimental work Wallace concluded that there 
were only minor differences in spectra computed with daily 
Or twice-daily data. Since the data examined in this study 
Ceo the daily variety, his conclusions, along with the 
tectetnat higher frequencies due to advection of short 
waves cannot exist in the model atmosphere leads the author 
to believe that the time interval chosen for this data re- 
cord was a sound one. 

Fn che BINED O2T program, the power spectrum is ob- 
enned by performing a harmonic analysis on the auto-co- 
warıance fumetıon of the quantity to be studied. Alfter 
Gaeh spectral estimate was multiplied by the interval of 
frequency (Af) of the band width, the spectral estimates 
Or Chat frequency band can be Summed to obtain the contri- 
eono thew tluctuations witmin that specific frequency 
band width to the variance of the time series which is 
under examination. 

Es US Otro pecan yestlega tion 

Thesp:rimarys point 8! interest in the spectral in- 

vestigation was point C (Figure 20) in the western bound- 


ary region. A thorough investigation of atmospheric and 
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Pecamte variables vas Conducted at this point with the fol- 
lowing results. 
a C daea at = point C. 

The data obtained from phase I (combined model) 
Cimthne experiment were examined first. The Fourier coeffi- 
cients from the atmospheric part of the combined model were 
analyzed and Figure 52 shows the power spectrum of A and B 
(solid line) and C and D (dashed line), obtained from Equa- 
muons 802.22) tom (2425), plotted on a linear scale against 
frequency (cycles per day). The power spectra of A and B 
Sree rdentical as are the spectra of C and D. The amplitude 
cite power spectra of C and D is slightly less than 104 
Mivthat Of Aland B. There is a pronounced peak at about 
yore days an the power spectra for these quantities. This 
Was expected since the phase speed of the imposed baroclinic 
disturbance in the combined model was estimated at 600 km 
per day, a phase speed that would yield a period of 8 days. 
In all likelihood, the estimate of the phase speed of the 
dats cur bance was slightly New. The 7.5 day period indicates 
that 640 km/day was the speed with which the disturbance 
moved across the 4800 km wide domain. 

Figure 53 shows the power spectrum of atmos~ 
pheric disturbance kinetic energy. There are some small 
peaks in the spectra in the frequency band from .08 to .24 
tesi oo days) but the major peaks are 
clearly at the lower frequencies. There are peaks at 


Ve aod on Ot eauOwwEoeeandmoe days indicating cycles in the 
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Figure 53. Power spectrum of atmospheric dist urbance 
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atmospheric energy on the order of one and two months. If 
seasonal variations had been included in the specification 
of the solar radiation and latent heat release, different 
frequencies of oscillation would probably be observed. The 
annual fluctuation in the above mentioned parameters would 
probably be reflected in the energy spectrum. 

Figure 54 shows the power pectrum of E, the zo- 
nal mean part of Yu thatetsmpeaoportional to the mean height 
and related to the mean geostrophic wind in the atmosphere. 
The quantity has a large spectral peak at a period of 32 
days and smaller peaks at 25 and 62 days. An examination 
of the cospectrum and phase between E and the disturbance 
energy shows that the disturbance energy rise at the 32 
faye period lags the increase in the spectrum of E at that 
period. This suggests that the mean flow increased first 
and because of barotropic amplification transferred energy 
Fomine disturbance. 

PERS er are per see power Spectrum of tire 
oceanic stream function. Peaks in the spectra are seen at 
periods ori 5 days (small), and 32 and 62 days (large). 
This is the first quantitative indication of an air-sea 
interaction phenomenon in the model. The 7.5 day peak in 
They power spectrum of the stream function is caused by the 
euer butLomwer the Awana 95 coelflicients to the disturbance 
component of the atmospheric wind stress (Figure 52). 

Since the disturbance components of the stress are small the 


Malos they cause are small and therefore the peak in 
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A omer spectrum Of ocean stream function at 
A pas s Period of spectral peaks indicated in 
days. 
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Mi pe tas. smaller at this frequency. The larger peaks 
atsabeutr 32 and 62 days are caused by the oscillations in 
the zonal component of the wind stress which is proportional 
tows. FAN examination of the cospectrum and phase between E 
amelssche ocean Stream function shows that the increase in the 
Dean atmospherie flow occurs) 3.5 days prior to the fluctua- 
Eaomssrusthe stream function at the 32 day period. Since E 
Rad spectral peaks at these frequencies it is not surprising 
that peaks in the same region of the spectrum appear in the 
ocean stream function., Figure 56 shows a 100 day time 
series of E and the oceanic stream function obtained by 
averaging the recorded data over seven day intervals. The 
averaging was done to remove the 7.5 day oscillation from 
the data to examine the longer period fluctuations in more 
detail. It is easy to see that the long period oscillations 
are present in both E and the stream function and that the 
maxima in E leads the maxima in the stream function by 
several days. 

Figure 5/7 shows the power spectrum of the zonal 


mean of the atmospheric surface temperature, Tye Since 7 
Ep roportLIomal to F, the peaks in the Spectrum are also 
nep esentacive ol peaks in the spectrum oit F., This was 
aer edib apectralily analyazime F separately, the plot of 


which is not presented. There are peaks in T, of almost 


2 
equal amplitude at 7.5 and 32 days, and smaller peaks at 
Dezronsyots23, 62 FanarIrlrdays.serisure 58 shows the power 


spectrum of Tea’ the sea surface temperature. Small peaks 
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can be seen at periods of about 25 and 32 days and a larger 
paa ia a period of about 62 days. These peaks coincide 
with some of the peaks in the Spectrum of pe Figure 59 
shows the power spectrum of the vertical heat flux between 
the atmosphere and the ocean. Peaks in the spectrum occur 
at 7.5, 32 and 62 days. The longer period peaks are at the 
Same period as the peaks in the atmosphere and ocean surface 
temperatures, indicating an exchange of heat due to a tem- 
perature difference between the two media. A comparison of 
the time series of the atmospheric surface temperature, the 
vertical heat flux and the sea surface temperature obtained 
by averaging over seven day periods is shown in Figure 60. 
It is readily observed that oscillations on the order of 
30 days are present in all three series and not so obvicus 
fluctuations on the order of 60 days are present in the sea 
ser trace temperature and the vertical heat flux. It can be 
seen that peaks in the sea surface temperature led peaks 
vertical heat flux which in turn led peaks in the 
atmospheric surface temperature. 

An examination of the cospectrum and phase be- 
tween E and I shows that the increase in E at about 32 
days precedes the fluctuation in sea surface temperature by 
Midas This was perno Sa SU creas ea tna 
mean zonal wind stress causing an increase in the wind 
apart the ocean circülation and an increase in the 
edvecwuonweh heat from the south at point C. An examination 


ormene wcospectrum of Ds EN and the ocean stream function 
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Sosa tre fluctuations in the Stream function at 32 
days occurs 7.6 days prior to the oscillation in T_ ,,- 
These phase relationships coupled with the previously stated 
eee eads to the following conclusions. The oscillations 
which were reflected as peaks Ta the pover spectra ‘Oh sthe 
variables were caused by an air-sea interaction phenomenon 
and were not merely spurious reflections of energy due to 
sampling or calculating problems. An explanation of the 
probable series of events causing these interactions is as 
follows. The increase in the mean atmospheric flow, re- 
cited ın E. caused an increase in the wind driven part 

Ene ocean Circulation 3.3 days Later. >This increased 
Smrcutation brought about an increase in the advection of 
mestzirom the South at point © and an oscallation in the 

Seis urtace temperature, atra period of 32 days, occurred 
das later. The total time from the increase in E to 
the increase in na was about 11 days. The increase in 
ea caused an increase in the vertical heat flux and an 
increase in the atmospheric surface temperature. An exam- 
nation trol the cospectrum andsphasexor See and T shows 
Enaertnesoseriliatıion in ee at) s2 days leads the rise in 
the spectrum of T, by 6.2 days. An exchange of energy be- 
tween DD and E is then required to maintain the cycle. 
Ensexaniwatcion of the cospectrum and phase of E and F shows 


Una e erpeak in F ©Oceceursapmior to,cherpeak in E at 32 days. 


When these interactions are taking place in the combined 
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model and oscillations at various frequencies are present, 
Fimosisevery dependent variable in the model will be affec- 
ted in some way. It has been shown that the long period 
eeni ations wiich are present in the coupled model include 
BiemoOcean as an important part in the oscillation cycle. 
Mou ld not be surprising to find that these long period 
oscillations do not exist when either the atmosphere or 
Secan is held constant in time. 

The data from phase II of the experiment were 
Seimed Next. With the time varying forcing functions 
removed, the ocean became steady in a short time with no 
Significant peaks in any of the spectra. Figure 61 shows 
the power spectrum for P Phi sirsi typical of all the 
ocean variables in this phase in that the power in the 
Spectrum is in the lowest frequencies. A comparison of 
Preure 61 and 58 provides a convincing argument that the 
eimesyaryine atmospheric eireulation is responsible for the 
transient responses of the ocean surface temperature at 
Ber loadssorsabout one and two months.: 

An interesting, and initially perplexing, pedak 
appeared in the ocean stream function at /.5 days in the 
acean only data (not shown). The ocean-only model required 
w enod Poiti me to adjust EO “the new forcine and during 
that time the stream function oscillated and then became 
constant. Although the energy in the spectrum was diluted 
Doe lone period ol constancy iie initial oscillation made 


Mee wevident at the frequency at which it fluctuated. 
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The data from phase III (atmosphere-only) were 
analyzed and Figure 62 shows the power spectra of A and B 
line) and © ands) (dashed line). The four quantities 
Mare ta pezk at a period of 8.4 days. This is different 
Eno the period at which there is a peak in the spectra of 
the same quantities in the combined model. For the 4800 km 
wavelength used in this study, a period of 8.4 days yields 
phase velocity of 570 km an This is about 12% slower 
than the phase speed for the atmospheric wave in the com- 
bined model. By eliminating the zonal variation in the sea 
surface temperature, the sensible heat exchange was changed. 
A comparison of the combined model and atmosphere-only model 
emows‘that the thermal wind was reduced in the latter. 
Bhıssreduetion in the thermal wind, shown in Table 4, was 
mieliwential in retarding the phase speed of the disturbance, 

The power spectrumser the disturbance kinetic 
energy is shown in Figure 63. A large peak at the 6.7 per= 
iod is observed along with a smaller, wider peak in the 45 
pomada ay period range. The peak at the 0.7 day period, 
which is near the 8.4 day period in the disturbance quanti- 
wee Ss perplexing in light of the fact that there Is no 
sharp peak around the 7.5 day peak in the disturbance kine- 
Phemenetcye(higuress3) in Ene scombinedemodel. A possible 
explanation is that the 6.7 day peak is very sharp and 
Vibe whe Nac a considerable amplitude, the total contribu- 


M oa he varianca oi the disturbance is not much greater 
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CATITUDE IHERMAL WIND THERMAL WIND 


(1) een 
12.0 0.0 0.0 
68.4 T 3.0 
64.8 Be) el 
61.2 og 35 
Dr Lt 4.3 
54.0 De 3.6 
50.4 6.6 7.5 
46.8 Ge) 9.4 
43.2 10.2 Se 
392 | Tia 
Sec 6.8 3.6 
el 4.4 I 
2808 30 SL 
IZ 2 28 
21.6 1 Do 
¡ESO 0.0 0.0 

TABLE 4. Comparison values of mean atmospheric thermal 


wind from combined model (I) and atmosphere-only model 
(III). (in m sec l). 
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Enemescher smail distuse indication of a peak around the 7.5 
dayeperitod in the combined case, In the eae Tor ihe 
quantities E and F, peaks at these same periods are ob- 
served, showing the exchange occurring between the disturb- 
ance and mean flow. 

Amsezanmınatıonsow data in phases I, II and III 
“=epotnt G indicates that theypeaks in the spectra noted ar 
about one and two months in the combined model depend on 
atmosphere-ocean interactions for their existence and are 
Hoti simply Spurious or caused by aliasing. The fact that 
these interactions are on the order of one and two months 
Penne che most significant point; because .the period 
would undoubtedly change with some other period of atmos- 
Meeris forcing. What is significant is that these longer 
podria oda clons can beMattributed to dir- -sea interac- 
tions. The absence of such peaks in the spectra in either 
the atmosphere-only or ocean-only phases shows that the 
fluctuations causing the peaks are phenomena occurring only 
when the freedom found in the coupled model is allowed. 
Although there is a peak in both atmospheric and oceanic 
variables in the coupled model at (7.5 "days, this is not 
considered as significant as the peaks at the longer per- 
Role because this Frequency vol oscillation was forced into 
the model with the selection of the 4800 km wavelength for 
Eresin le mave inex. JGhancingwihis mraveleneth would cause 
Micmpedkewingethe speetrim at 7.5 days to shift and undoubted- 


ly would have many other effects on the model. It would 
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ere stine Eos nene peaks simsthe Spectra at the 
Tower frequencies would be shifted. 

Spectral analyses of the kinetic energy in the 
ocean of the combined model at the three ocean levels were 
meat point C. An approximation to the total kinetic 
Meme edt point © was used for this study The total kine- 
weee nergy should include the vertical mean kinetic energy 
and the vertical shear kinetic energy. A mathematical for- 


makarıon for this is 


KE ži! y A ÓN Sty?) (6.2) 


ocean 
In recording the data at point C no provisions were made 

i ow 
tonr the calculation of EN and so the zonal component of 
Moe vertical mean current at this point was neglected. The 
Zonal component of the vertical mean current was small and 


Buererore this was not a bad “assumption, ss the kinetic 


energy computed at each oceanic level was 


KE 2 st (ut)? Paes 2)°], (6.3) 


ocean 

Figure 64 shows the power spectrum of kinetic 
energy at 100 meters in the ocean of the combined model 
Versus sirequeney in cyceleszper dasssborkeplottea on Log 
scales. Figure 65 is the same plot for 500 meters and 
Figure 66 shows the 1300 meter depth analysis. Peaks in 
Co ect raw atsperiods Of 07. 25-32 eames.) days are. ob 
served at all three levels. This is consistent with earl- 


ier findings. The amplitude of the 25-32 day period 
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Figure 64. Power spectrum of oceanic kinetic energy at 100 
jot pot Cn phase ls. Periods of spectral peaks 
imdreated in days. 
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Figure 65. Power spectrum of ocean kinetic emeroyvs27.500 
Poo wee point C in@phase I. Frersods of spectral peaks 
indicated in days. 
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che largest at the upper level while the amplitude of 
Ene 62 day period and 25-32 day period are about the same 
at the two lower levels. A comparison of Figures 65 and 66 
shows that the spectra for the two lowest levels in the 
model nearly coincide. The maximum amplitude in the power 
spectra decreased by an order of magnitude between the 100 
meter level and the 500 meter level. This lessening in the 
variance of the kinetic energy between these levels is not 
Surprising since the lower levels are less active. The 
amplitude remained about the same between the 500 and 1300 
meter levels. 

Thompson (1971) analyzed a long series of data 
recorded by the Woods Hole Oceanographic Institution at a 
site esca at 39 200R OW in the Nori h Atlantic. In the 
model this would be analogous to a position approximately 
200rkm northeast of point C. The data were recorded at 
900 second intervals. The folding frequency associated 
with this small recording interval allowed Thompson to in- 
vestigate the shorter period oscillations. Some gaps in 
the data forced the estimation of some of the spectra in 
Preto tre queney. region He analyzedithe northward and 
eastward components of velocity and through a summing and 
averaging process calculated the "kinetic energy density" 
at several depths in the ocean (not shown). None of these 
spectra showed any peaks at periods greater than one day, 
although osciallations with periods of about ten hours were 


observed. The spectra for the three lowest depths were 
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edt coincident, and there was a large decrease in the 
amplitude of the spectra between 100 meters and 500 meters. 
In this feature the study of Thompson and the study con- 
ducted in this work were similar. 

DS Cua y —Oted2tanateypoint A. 

At point A (northeast area) the power spectrum 
of the coefficients A,B,C and D were greatly reduced la 
maximum value of 16 CA cycles per day for A and B 
aesechmpared to 675.0 at point C], although the peak in the 
Spectra was at the same period, 7.5 days. An examination 
orsthe coefficient, E, shows the spectral peak at 32 and 62 
days and a large peak at 120 days. Looking back at Figure 
33, a Smaller peak in the region of the 120 day period can 
be seen also. The oceanic stream function at point A shows 
oscillations that have smaller amplitudes in the power 
Spectrum than do those of the same quantity at point C. 

The biggest difference between the power spectra of vari- 
ables atipoint A and those =atepointec occurred with an 
Fete roe nO peaks at thes 201, O2edayeperiods and all the 
power is in the lowest frequency, much the same as seen in 
Figure 58, the same quantity in the ocean-only case. 

AC point A in theroccan only phase, the power 
Spectudmeccemb led those from the saner phase at point C. 

The maximum of the power spectrum that existed in the low- 
eSterregiiene, Yvas increased put villinot be discussed fur- 


ther since it adds no new information. IN Eee aso Ma Ene 
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Doreen theslowest frequeney merely imdicates that the long 
term change was larger. In the atmosphere-only phasc, the 
ES peca cor the Fourier coefficlients/’resemble those 
Ac except that the maximummor the Spectra at the 
8.4 day period is reduced to approximately 10% of the value 
o ed at point €. At point A, the meridional temperature 
gradient is rather weak and the amplitude of the baroclinic 
wave is much less than in the lower latitudes of point C. 

coo cudy of datagatecoint Bl) and E 

At point B (central area) the power spectra of 
emesecoefficients A, B, €, andi D were slightly larger than 
eneee ral point C. This location at 43N is in an area of 
large thermal wind, so the amplitude should be greater. 
Except for the increase in amplitudes the spectra are very 
similar to those of the variables at point C (Figures 52 
through 56). The power spectrum of Tea’ though, is more 
like that at point A, with all the power at the lowest fre- 
guwemey. nis indicates extremely Tone period fluctuations, 
Po pablito long for this data record to portray. 

The ocean-only variables at point B resemble 
those at point A and will not be discussed. The Fourier 
evetiacitents, A through DG Nave a pealyat the same 8.4 day 
period and the maximum in the power spectrum is increased. 
No new information was obtained from point B. 

DSeparateldiccuss iongo ea ta at point D 
u no e gmade because ol Soini laney results to point 


eo oe aee eul ts were obtained when the data at point E 
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were spectrally analyzed. dhe Fourier coefficients, A 
Ao Dad multiple pesks an the speetra that were not 
representative. It was decided that this point was selec- 
BodmEoo close to the southern boundary and the artificial- 
ity of the wall caused unusual oscillations in the 
atmospheric data. The results from point E were discarded 
waith no further comment. 

The results of the spectral study were that the 
ocean circulation caused transients to appear in the atmos- 
phereic circulation in the one and two month periods and 
possibly with a period of 120 days and that similar tran- 
sients appeared in the ocean circulation due to the pre- 
sence of the fluctuating atmosphere, A great deal of 
information has been wrested from the spectral statistics 
discussed here, but a great deal more could be obtained if 


time permitted. 
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WERNE SVDNARY 


A numerical model of a combined atmospheric-ocean system 
masmdevelopbedsior use in tan Investigation of the long-term, 
large-scale air-sea interactions and the effect of these in- 
Benacrjonssonsrransıemteszthat appear in the atmosphere and 
Oceanic circulations. The problems associated with observa- 
tional studies of the air-sea system were eliminated by 
acios 42 numerical, model thal could produce a long-term simu- 
lation of an interacting system in a reasonable amount of 
computing time. 

The atmospheric model used in this study was fashioned 
after the two level, quasi-geostrophic model developed by 
Phillips (1956). New dependent variables were introduced 
by the use of a one-dimensional (in x) Fourier series, in 
which a single mode was retained, to represent the origi- 
wo variable T lnis method of approach allowed the rapid 
integration of the atmospheric equations since the number 
of degrees of freedom were greatly reduced. Without this 
modification the combined model could not have been inte- 
ei telsTorgarsusricesenely Long time period to obtain mean- 
ingful results. 

eno cean model utilized in this work was a modified 
version of Haney's (1971la) primitive equation model. The 
modificat iona included a reduction in the number of levels, 
a oran tal region size (and the number of 


Atmos mae Teduetton in the depti of the basin and the 
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conversion to ß-plane (cartesian) geometry. The primary 
reason for these modifications was to increase the ratio of 
simulation time to computing time and permit a long term 
simulation which was accomplished. 

The coupling of the model was accomplished through the 
exchanee (Of momentum through the wind stress and the ex- 
Change ol heat through the turbulent exchange process which 
was proportional to the air-sea temperature difference. 
Assumptions were made to simplify both of these exchanges. 
ie atmospheric model, however, had to be restricted to 
Simplified formulations if it was to be integrated success- 
fully for the necessary length of time. The absence of 
Beristeure in the atmospheric model required further simpli- 
rear cn ein the formulation of the latent heat flux. These 
Simplifications were indeed necessary and an examination of 
the results showed that the model produced atmospheric and 
Oceanic circulations that were adequate for this experiment 
in their portrayal of the real atmosphere-ocean system. 

The experiment was conducted in three phases. Phase I 
was a long term integration of the combined model. The 
computing time was reduced by the use of a larger grid size 
orale irst portion of thrszpnases Tine 400” km grid dis- 
tance was that used to accelerate the integration of the 
combined model equations in the first part of Phase I was 
Just adequate to resolve the western boundary current. 

Picco neient ofr lacera dillusi on of haat had been re~ 


ee o reven t large horizontal transports of heat by - 
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sion tanda Keep the meridional temperature gradient at a 
as ticiy Men level. The lateral@diffusion of momentum 
was reduced, but could not be made small enough so that a 
Guu on ideo Lotion eould be obtained without having 
hendt ao eaen boundary current become smaller than the 
we onta erid distance SE he unrealistic spatial oseilla- 
tions seen in the results obtained with the 400 km grid were 
directe result ol the Farge zrid Size. 

The spatial oscillations observed in the 400 km grid 
model results quickly disappeared in the 200 km grid version 
of the combined model. A state of quasi-equilibrium was 
reached in the long term integrations of the combined model. 
The results shown throughout this thesis were considered to 
WMenSatiStfactory for the purpose of investigating the effect 
of the oceanic (atmospheric) circulation on the time-depend- 
ent behavior of the atmosphere (ocean). Some features in 
the results considered worthy of note are summarized below. 

A western boundary current developed with very strong 
northward flow, analogous to the Gulf Stream, on the west 
Soto cie sub=tropical gyre and strong southward flow, 
analogous to the Labrador Current, on the west side of the 
sub-arctic gyre. The gyres themselves are positioned as 
found in nature and are a result of the surface wind stress 
developed in the atmospheric model. The conditions in the 
HO@EleadstecoOuner of the domain are similar to those found 
ana nesnortheast Pacıtic. (new warmeretcater is due to the 


SoUvemwesence Of Sane vertical shearfeurrent into the 
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easter bora resulting in downwelling in that re- 
gion. The region of warm water and weak currents in the 
SO ent ral recom oi the domain is analogous to the water 
Faumd un the Sargasso sea area of the Atlantic. 

The zonal winds in the atmosphere develop a horizontal 
IM taco iS cs e similar to that found in the real at- 
mosphere. A three cell meridional circulation was formed 
and slowly meandered several degrees of latitude north and 
south of its mean position. 

The second phase of the experiment consisted of a six 
year integration of the ocean-only model and a state of 
equilibrium was quickly reached. The upper level of the 
ocean had reached a state of almost steady temperature but 
the lower level still showed a very slow increase in the 
temperature, This was a good example of the extremely long 
time necessary for the vertical diffusion mechanism to bring 
the deep water into thermal equilibrium. 

Phase III of the experiment consisted of the integration 
of the atmosphere-only model equations. Energy balance 
studies of the atmosphere from this phase were made and com- 
pared with the same studies from phase I. The energy 
studies compared favorably with each other and with the 
existing knowledge of atmospheric energy exchanges. 

Datawproduced by the combined model and recorded at 
S ee Mounts weres spectrally analyzed yielding peaks in 
Bu ze er ao Bothrocesnsic ind atmospheric variables at 


periods on the order of one and two months. Cross-spectrum 
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anaty ses and time series analyses of Several atmospheric and 
oceanic variables led to the following conclusion about the 
air-sea interactions. The longer period oscillations were 
Pre entimerhe data andiye attributed to air-sea interac- 
tion. An increase in the mean atmospheric flow preceded a 
mse in the barotropic flew in’the ocean by about 3.5 days. 
This increased oceanic flow caused an increase in the advec- 
tion of heat and oscillations in the sea surface temperature 
about /.6 days later. This change in sea surface tempera- 
ture led an increase in the vertical heat flux between atmos- 
phere and ocean and an adjustment of zonal mean atmospheric 
surface temperature. It is interesting to note that even 
with the limited model developed here, there are longer per- 
iod interactions between the atmosphere and ocean. If the 
model was modified so that more freedom were allowed, then 
wmeeractions at different periods might occur. There are 
several ways in which this model might be changed for future 
Studies. Some suggestions are listed below. 

The incorporation of continental land masses and the 
corresponding allowances for a zonal standing eddy into the 
ode szacnecessiıty. Miis erfect would modify the heating 
im the wmedel in that 2 morezrealistic longitudinal pattern 
of vertical heat flux would emerge. The effect of this 
change on the frequency of response in the model would be 
Tite esoo seres neo pp oration of a second har- 
monic into the atmospheric Fourier formulation would simu- 


mate thezeontinental effect. Inis vould also cause a shifr 
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in the position of the short period oscillation or adda 
new period of interaction. A study of an atmospheric mode]. 
with a two-wave structure was carried out by Trumbower 
Gaye Seciowineg the feasibility of including a second wave 
w e abmoscpmerlc formulation. ineorporation of the second 
Wave in the atmosphere would require an increase in comput- 
nee immer Dur seomparedstostheztotar time required for the 
experiment, it would be negligible. 

A change in the sensible heat term should be incorpora- 


tcd in the model. "in the term Ha E. - T,)5 the quantity 


ea 


Hy was held constant. This quantity should be divided into 
wo parts, one that is constant and one that is a function 
Simcoe surface wind speed, This would cause a portion of 
the vertical heat flux to be dependent on the wind speed 
and allow more freedom in the model resulting in other 
Worms of Interaction. 

An analysis of the energy budget of the ocean should be 
welor porated in any future calculations. The transforma- 
tion of energy between the vertical shear kinetic energy 
audstcthezswertical mean kineticwenergy and the dissipation of 
energy is an area well worth studying. 

themoiases TII experiment should bel carried out again 
a aur ace temperature with an x dependence inclu— 
ded. A zonal variation in the heating would be present 
and the frequency at which the yariables oscillated may be 


changed. 
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Domzıcenoephersesprediertonstor periods shorter than a 
few weeks, an interacting model may not be necessary. An 
atmosphere~only model with an accurately-observed air-sea 
temperature difference may be sufficient. However, longer 
Yange predictions on the monthly or seasonal time scale 
Should include some kind of simultaneous ocean calculation. 

The atmosphere-ocean model was synchronously integrated 
Does year si Mulat lonor he eireulations of the two 
media. The interaction observed at periods of one and two 
months provides evidence of the importance of the effect of 
the ocean circulation on the transient behavior of the at- 
mospheric circulation, and vice versa. Other periods of 
interaetıon surely exist but the investigation of all of 
them is beyond the scope of this work. It is hoped that 
the results presented in this paper are sufficient to moti- 
vate further research in the important realm of large-scale 


air-sea interaction. 
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RERENDIX X 


Derivatıonsaresrauerrons. (2,22) to (2.27) 


Bauart home (2 nlc ie Mayan ibewmiew rl fen 


o „2 a ay o 
EY Wy? + Ik x VW ag V CV Yay? aller vv. 
, ay oT one 
V(V vn) T B Vx rare A A ~ dy in , (Al) 
where 
UY E co x. By, t) Sin kx, (A2) 
Ver IE AO y A coses F Dy, t) sin kx, (A3) 
GE O W, | u, (A4) 
and 
er Zn w, | Vy: (A5) 


Ay examination Of quantities. in Equations (A4) and (A5) 


shows 
OW), 
E = dy , 
ov, 
Me eae ® 
1/2 
ow, 2 ow, 2 
and lv, | [Ge ee) : (A6) 


linear interpolation wadas used to obtain 
3 I 
Py = o Va => Va , (A7) 


and the relationships 


E 
Vy = > > 

A 
Vy = a 
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and 

Ya? Uy = Vo: (A8) 
Then 

Y, = Vy 20 - (A9) 


A Equ ions CAS anda AI) Eo (Ab) to obtain (wien 


2k and higher order terms are neglected) 





__ 28, 4 aF_ (aA _ » ae 
u, = Sy + 2 Sy 57 2 EZ cos kx 
OB 9D l 
-- T - 2 Jy’ sin kx, (A10) 
kB- 2D) cos ks ek 7 Ke): sinnkx, (ANI) 
and 
meea | 1,94 _ ¿20,2 , (dB _ 52D)2 
WV | | ay - 25) | m 253) + (ay ee) 
Cec (eee | 
JE OA EE JA OF 
ty 3y 7 “oy By 7 “ay dy 
+ 800 ar 3E 3B _ ,dE ƏD 
: n T cos kx F O y Po T 
t72 
OF OB ID OF j i 
= Er Er T I dy? SIn aa : (A12) 


Equation (A12) is simplified by making the assumption 
that contributions of the disturbances are much less than the 
contribution of the zonal wind to the magnitude of the sur- 
face wind. The cosine and sine terms in (A12) are therefore 


nepieetedsand a new formulation for NV, | istobtained: 
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JE E > 
AEE 2 Ge Za) 
1/2 
2 
+ > E an Den) | o (A13) 


Equation (A4) is now 
] a BE , „BE, _ aA _ z2C 
ce PC, NW, | fr- — + 2] y Be GOS) kx 


OB oD ; 
-- ae - nn sin kx > (Al4) 


and Equation (A5) is 


a pc, AV, | [am cos kx - k(A-2C) sin “| ; (A15) 
oT OT, 
The term = dy l4 can be written 
OT oT 
Be Xie LA oF _ 2E 
lax 7 By l4 7 PSp ag ay Sy! 


+ LA k“ (2C-A) - AN C - m | cos kx 
+ A k2(2D-B) - O 2 = | en xf aa) 


When Equations (42) (A3) and (ATO) are substituted in 
Equation (Al) and terms of wave number 2k and higher are 
neglected, an equation of the form 

—[G E) y Es (Y E oO G ayt) sin kx] 

a 2 3 

= F, (y,t) + F,(x,y,t) cos kx + F¿(x,y,t)sin kx (A7) 
is obtained., Coefficients of the terms independent of x, 
coefficients of the cosine terms, and coefficients of the 


sine terms are equated to yield the following three 


equations; 


PC 








we sen) ger, 





Q 
—| Kae ee ee 
ot Jy? dy 9 y? O y dy? dy 3 
-3 y, > : 
Ea eg By DY] A Bk 
dy 
2a Y ne oa ies oc _ JA 
+ = scp] IW, | k“ (A-2C)+ ay l Vl CS ay) 3 ‚ (Al8) 
2 2 2 3 y 
ui E 
dy? y Jy? y dy? dy? 
ee (Caer Ci] ee Ak 
ae. V E Ce „eD _ OB 
i Ar ey | Ny k“(B-2D) + ay EI, | C255 jy)! , (A19) 
and 
2 2 2 2 2 
Es: E 22 ee O moe) 
e le es ee a 
4 5 0, MIO = 5.) 1: (A20) 


Rewrite Equation (2,15) to obtain 


o 
ee le yo), VO, 





oy 2 
tik x Wy e VCV? +B yo en 
O 
ont OT on OT 
A 2I El - [2 - 1,5, (A21) 
Ap ðX yee OX dy 74 
and examine the term HERQ/2E,. Substitute for Q to obtain 
: 2f 
w*kQ _ pkg o oT 
z OT sen orem et COGN Ge) CA Zz) | - (A22) 


NIY 





Define a constant, C} such. that 


u IT. R 
ee =e emer Ca A Dn (823) 


and the term is then 


2 2f 
Paca SO: 
26 8 p eee RV + 017. 2 


Substitute Equation (A3) for vn and 


DIE 
2. 2 
YKQ _ UK ee 
2£ | a: p Hy en R [(F+C,)+C cos kx + D sin kx] 
(425) 


is obtained. 
OT, OT 
The we rio) np 2157 - ay 12 is then examined. Recall 


thart 
D RE (A26) 


totr can belbroken up into its vector components 





z Ps = 
"x, (pk), ( or (u, uy) (A27) 
and 
2 ps E 
Es (pK), ( ur (va = v). (A28) 
The quantity T, Pay be written as 
OW JW 
2 pe A a 
aor (pr), ( D2 ( m, (A29) 
that leads to 
dT d Y 
ee E (A30) 
oy 2 p eo 
y 
The quantity To may be written 
ou OW 
= pe ee. 
O A A (A31) 
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which leads to 











IT dY 

PSA PE I 

37 2 (PK), ¢ ay a (A32) 

x 

and 

oT oT 

_ X aa . PE 2 

Finally 

Br en Be On ee 

PAS 2155 Tae Lp ies 4 (ok), ( on V Ver (A34) 


and can be written in its Fourier components as 


2 2 3 2 
Mo = (Ck -= a eos e 9 D 
dy 


gy oy 








DiS ne 


-6 
Ap 
(A35) 


Repeat the procedure previously used on Equation (Al) 


on Equation (A21) to obtain three equations. These are 
226 2 
15 - 8°C -n?c]= k 1 97D y CE oe SE yc? + p?) 
3 
3 SF Bat - uf) - os a BJ] - PB, Dk 
y 3y” gy? O 
AEE a g PE DE 2 
O d y 
Bor = oc _ oA 
+ A p TOS -A) o, vi ; (A36) 
9 ¡9D | a 
el, -MD-u?DJeK[sock?+u2)+ SEA? - u?) 
2 dy 
d y 
IE 2 2 n 2 
y dy? Y dy? dy? ay? 
°K > g pg 9D 2 
a A RON TS kD) 
er V k? IDe BIE DT 9D _ JE 


Z0 





2 2 2 
2 - rk) 2720, _ „ee 
y dy? dy? 





+ IB 2 DIA a ae an, 
Jy? 9 y? dy 
> De 
P 2f Ê p e ae FR E TE 
Mot, Els e en. (438) 
Ap 2 p 2? Jy? Ap Dey A oy 


This completes the derivation of the equations for pre- 


dicting new values of the Fourier coefficients. 


202 





ARBENDIX DB 


Derivation of wand Vy Bauat3ons. 


Applied at level two, Equation (2.12) is written as 
W y ae 
d 3 1 3 1 o Be 


Sobereitute from Equattonm(2 io which yields 





ð : oap MEE. € 
zt Pr tik x VOY Vin 2f Yo a Qo: (B2) 


Solve for U, which gives 


2 


O A ae 
2 oA p ot oy ax Ocean: 2f 2 i 








Subst urneerrauer one ale) andy (2.17% into (B3) and 


neglect terms of wave number 2k and higher which gives 


27: 
Be cr a on 
W, = o Dee) + u a 
KQ 
OF 2 oc JE OF 
+ a DE ele - Aya - Ol cos kx 
oD JE OF : 
ie T S — jy") J sin kx 
or 
0, = w+ (0,) cos kx + (w,) sin kx. (B4) 
Tew wan tay Wy is composed of a mean component and a 
disturbance component. The mean component may be written 
-fo IC DA 9B aD 
Du a 02 > Sale A 
a, non en 
KQ 
OF 2 
es a ey 
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aıdscrwe dAisturbance part written as two components 


ou 
u O ðC E, r 
A 2 
and 
Zi 
_ O dD Cenacle 
ea ya eG) (B7) 


The mean meridional motion component can be derived 
irom the zonally averaged continuity equation. This is 


writ en 


Vy oe = 0. (B8) 


= (B9) 


and this is solved numerically to obtain Va as any polne in 


MC MOde! walter obtaining W 9 from (BOIN 
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AERENDIX € 


Beriyatıion ol Energy dand Energy Transformation Equations, 


Pomudetive cherenereysrelatienslips, the A, B, C, D, E 
mdf wequattoncm|(2.22)eto"ce,2/)) were multiplied by A 
through F, respectively. After further mathematical mani- 
pulation, the energy equations were obtained. 


Multiply the A equation by A to obtain 














2 2 2 3 
ass A ei E ka [$E O° 3B z — Jen En 
y dy? y dy? dy? 
3 
A o - Bo aBk 
A 3 ys dy O 
BA 
E. 9 oc _ 3A 
+ E pc, A WE ee 2 C ERE NAO a a) | (C1) 
MAN Re B eguation by Bo and add to Equation (CCl) to get 
2 2 2 
ASS -o ake BÈS -o Bk? = kpa 23 
dy" dy? en 
- 2 q2 3 
E poe 97 aA A ask 3°D _ oF. aC _ ADS F 
DÁ dy 2 y dy? y dy” dy? 
3 
+ Bee - eo eau DA - au CB) 
Jy? Oy ay 
= 2 2 2 
A > pc, NW y | k [(A’-2AC) + (B -2BC) ] 
+ ASSN, L5G - 3] + 35501 e re (02) 
4 Oy 
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Treat the C and D equations similarly and add to obtain 





, 2 
1 - ck” = cur] + Do „22 - DI ED bs 
t 2 
dy ay’ 
2 2 2 2 
o - pa Cy Me ct bu pa A, 
es dy Zar ay? 
3 
- Lac - a0) - (x? u?) PE (BC - AD) 
dy? 
— 2 
+ 75 PC, Im, ERST G2) oe -c2 lv, 10255 = Er 
9 = 9D o B g pg dee 2 
Ep AÑ 4. 4 aa Ck 
DyN, ¢ > u Er (PK). ¢ pale, ) 
2 
ee een (03) 
dy? 


Coubine aquattons (02) andecG3)eand integrate over the 


Her on tO obtain 











y 
2 ‘ 2 2 
J AZ E BL o] 0 A 
dy? dy? dy? 
O 
D 2 2 A JE, 3 B 9% A 9%D 
Deal ak -Du ]\ dy- nf le Si = 
dy? yay? by? dy” 
O 
2 
poe 
dy” 
2 2 2 2 
+ nl Do po A BI +0) 41? (Bo - AD)] 
2 ay? dy? dy? dy? 7 
Ron _ TE oe 2 990 _ BA 
a pc, {CA ©) (7, 1 (A-20)+ Sob |W, | (2 nn 
ta-») (7, RA TER = >) 
g pe 2 9°D 2 
+ N 4(pk), ( = „te -Ck )+D(—— : DE) ] > dy. (04) 
dy * dy 
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Perform the multiplications and apply the boundary con- 


dirtons. to acquire 


oe 
ot 


f ds, 24 (BB) 24 (BS) 24 (By 24 02 (A245 2407402) 4u2 (C242) ay 


2 
po C 
y 


2 
g 


ae 


dy? 


2 2 2 2 
I P a P De 
O y 





2 2 
y BF (,9°D _ pala 


Ay Jy? Jy? oy Jy? 9 y? 


BE — | 2 
2 y (BC- -AD) ]- m PC, [cae (li, 1x (A-2C) 


ay ENV, Cage + 5791 +0) (NW, [9 (3-2D0)+ 511%, 1 (257 


O (= EL a y ( 22D -Dk2)] 


dy” dy” ja 


(C5) 


Apply similar techniques to the E and F equations to 





Obtain 
y y 
2 2 
3 = 2 at pee ae on oD 
E zy ) + +u“F ]dy = ie B E 
O O 
2 2 2 2 
m a bp Be- -AD)] 
gy? Y ay? ay? By? dy 
2 ZE 3 = 
= TEN T) S y 
E Zu Een p pl ¢ Er | 4 | 
oak ro A ay. (C6) 
9 Oy 22 ay? 
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Le zezomalsipokenkial energy, 
Po PECdayepovemrela l ener y , 
K, = Zonal kKinetierenerey, 
KR = eddy kinetic energy, 
€ = conversion of P O 
p Z e 
C = conversion of P COCK 
e e e 
Cr =z Conversion ol K EEO; 
e 2 
E = Conversion oi KO P, 
Z 2 Z 
G, = generation of Er by 
D pee dissipation o RD 
ei e 
Dos = dissipation of Ko by 
D = dissipation of K, by 
and 
D = dissipation of K_ by 
ZS Z 


drabatic processes; 
internal friction; 
surface frictrion, 


ineernale<efrıctien 


Surface friction. 


From Fhillips (1956) and the basic definitions Tor Ve 


and UY 


used in this work (Equation 2.13) equations for the 


potential and kinetic energy components can be written as 


dy = 


fo 
fi 
fr. 
/ 


dy = 


75 


dy = 


and 


2 
2 
BR dy, 


(C7) 
Eje + D^) dy, (C8) 
a 
ON SP Jay, (C9) 
[ich a? ES + 
Peele? 
ne (A +B °+C°+D aa (C10) 
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The right hand sides of Equations (C5) and (C6) yield ex- 
Pressyons for energy generation, Conversion and dissipation 


as follows: 


2f 
a Me _ O 
G, | re PE A u eae bey, (C11) 
Be Deo u cay 
t = | [> uo CBC- AD) dy! dy, (C12) 
: 2 2 2 2 
= - leas E A 
\ | y dy? dy? dy? dy? 
2 2 2 2 
a ne Ze lay, (C13) 
ZN ay? ay? dy? 
2 2 
ee e S T 
ei Ap 22 Jy? 9 y ? 
er. n P Ie a er oC _ 9A 
De Fla PC, K o) (1%, E : Iw, | 35 a) 
+ (B- o> (I, x? (B=2D)+ 51,105 > = RECS) 
= — ie pe ax 
D Eo oa UF T (C16) 
y 
and 
2 ces a ON A 
DE Lis pC, E F) aylly, Ge >]! dy. (C17) 


The quasi-geostrophie first law of thermodynamics, 


Equation (2.9), can be manipulated in a similar manner to 





deriye expressions for Co and > These are 
2£, = 
Fer d > ¡A (C18) 
and 
oe 
Co = J IT, (wc Ar WD) J dy. (C19) 
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NN DES 


PRcwdmeeTrerence Lorm ofstheTeauatıon for the prediction 


of the vertical shear current, u', analogous to Equation 


es 
zn» - oo Q ie Q 
Pee Kg we 
2At Ae ane 5 an nr 
K=1 
(L+1) (L-1) 
' = i 
+ a TR +(1 Sn ) (D1) 
where 
L _ I _ ees 

isk o (Az), Sei 4 Mig Ma jr 


(u) L = (w, „tw 


13 913 Mig, gd k+ Mz | 





1J 
i 2 
A A a pg, 
+ —— |(d2), A [lo E u a) 
(Az), k “m Wigig 7 Mag? ag 7 Mim] 

A A a 
Al dis E er oe 

(hz), 1 (Az), 3 


and 
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and 


AA = Ax Ay 
and 
q A UK a en) 
ee De ee 
re je 
= Ea + y ) an (u + u ) 
IL ij 8 ij ij-1 
A 
= any kl 
1 Ax (Az), 
ee C D ee 


A ze O a y) 


1 Ax(Az), 
A i ee (ov, Pacis? 
a Or) 


(D3) 
The difference form of the equation of the meridional 
1 


component of the vertical shear current, v', analogous to 


Equation (3709) Ls 


0 ame Ce 


e E km 
ik Tjk oN A > £ 
2At ee, nz SE (Az), 
k=1 


2E 





(241) Cam 
p3 1 1 
E a ] (D4) 
where 
A! AA 
Sn (Az), O a Ud kt 


iz Os Marg ages Mata g tL ety et aj 


g 
A 


1 
Ax(Az), GE ng O = ij 5191 


1 
Jä 


O 


1 
= Gar A a 


+ ( - (v 


O UN 


N ea See e Kae] 
(Az). 1 Che eee 
=a => 


a 


k(v ye k (v eye 
| (D5) 


l SE 
2igk “pleas agape’ TE (Sir; a 


+ (u. 


o Dr + (u 


ee er: Zn, 
Ay (Az), 


1 
AN a5” lt E NETOS 


+ ( + u) + (u, 


E 0419) 

Ax (Az), 
eh. ee I 
et 8 Ic rl 


“i3+1 


-{- 


[5 (v ) 


Ze A Se ZES 
AX (AZ) 


1 k 
NA is? ee Vi-13) 


Y . (D6) 





Or 


202 





At each time level &, Equations (DI) and (D4) are solved 


Cori 


tagether for the unknown shear currents, ene and 


asi ze 
ijk 
The difference form of the vorticity equation analogous 


to A C3. 1e) is 


/ 
LCP) by a TARA HD )- 


1 1 

Z A o s PEN. 
l 
2 





BGE) RE 
E itl) A 
t Por ag Ea )] ou 
where 
Q eee i | Q 
a, a i ED age #(A2), AylCagyy 58g 5) 
k=1 
Coe) 
A A er 
k(u zu) k(u - u ):. (2+1) 
k-1 "Ki ko "k+1/i 
+ AA [ 2] (D8) 
(AZ) Giz eo pa i 
k-7 k+> 
km 
2 _ 1 1 Q 
San ti TO o A 
a 
Cosy) 
ee Den este 
k(v - v, ).. Sr jr 
en. 8 ij 
“2 
and 
= Ł 
A (ae 


Equation (D7) is solved by relaxation techniques to 


obtain T. 
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The difference form of the thermodynamic equation analo- 


cousto Equation (3.16) is 


(2+1) (2-1) 
fae er ) 


2At A N 


(T 
ijk Oe ES 


= +w.. 
“Yij+l Ť “Yi+1j+1?k-Ż Du a 


2 
E Yi+1j+1? k+ž CD ket A E eee) 


2 l (2-1) 
A ee, 


Typen Tu? EO Tee?) (2-1) 


De oc En cc le 
ey, 1 (Az), ij 
2 a 


+ 


(Az), 
eet (ST). (D11) 


-- 


where 


IH 


ren 


ray 
< 


A 


and 


(Az) 
nn! (ST). is the convective adjustment discussed pre- 


viously., 


Zee 
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